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ABSTRACT 

Chiral restoration phase transition in hot and/or dense hadronic matter is dis- 
cussed in terms of the BR scaling based on chiral symmetry and scale anomaly 
of QCD. The precise connection between the scalar field that figures in the trace 
anomaly and the sigma field that figures in the linear a model is established. It 
is suggested that in hot and/or dense medium, the nonlinear a model linearizes 
with the help of a dilaton to a linear a model with medium-renormalized con- 
stants. The relevance of Georgi's vector symmetry and/or Weinberg's "mended 
symmetry" in chiral restoration is pointed out. Some striking consequences 
for relativistic heavy-ion collisions and dense matter in compact stars following 
stellar collapse are discussed. 
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1 Introduction 




One of the most intriguing problems in physics is how chiral symmetry is restored in 
matter as it becomes hot and/or dense. There is presently a lot of controversy and there 
are many interesting suggestions. 

Some time ago we proposed that scaling, which is a property of the Yang-Mills 
equations at tree level, applied to the dynamically generated masses of hadrons made up 
out of the chiral (up and down) quarks. Our result for vector (V) and scalar (a) mesons in 
the low-energy sector was 

11 ~H$L S3 ••• (1) 

f n ray m a 
The nucleon effective mass scaled somewhat differently [§] as 

(2) 

In these equations the asterisk stands for in-medium quantity. 

Since our results were obtained by introducing the breaking of scale invariance with 
the low-energy chiral Lagrangian through a scalar field denoted Xi which we called the 
"glueball" field, our paper has been largely interpreted erroneously as tying the scaling (jl|) 
and (|2|) to that of the gluon condensate. We wish to emphasize that this is not so and 
clarify where the error in interpretation is made. In our argument in Q, eq.(13), we split 
the glueball field into 

X = X* + X' (3) 

where x* is the smooth gluon mean field and x' is t ne fluctuating scalar glueball field. This 
separation effectively splits x i n t° the field x* which scales the quark condensate and hadron 
masses, and the "non-smooth field" x' which governs the scaling of the gluon condensate. 
In the role of the latter was not clearly specified. The more precise way in which the 
separation, guided by some results of lattice calculations, is made was given by Adami 
and Brown ||] who brought out clearly that there are effectively two scales in the broken 
symmetry sector of QCD These scales can be given in terms of the two "bag constants": 



1. The bag constant for chiral symmetry restoration 

B x sb « (140MeV) 4 . (4) 



* 2 Several other authors have arrived at different conclusions. See, e.g., |E| and footnote #12. To 
distinguish the scaling we are advocating from those favored by others, we will call ours "BR scaling" in this 
paper. 
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This bag constant is essentially the condensation energy of the vacuum in the Nambu- 
Jona-Lasinio (NJL) model ||. Basically, in this model this is the amount that the 
condensate of negative energy constituent quarks is lowered by developing masses. 
This discussion is, however, model-dependent and, while the physics seems sensible, 
it cannot be satisfactory from a formal point of view. Our aim is to improve on it. 

2. The bag constant for gluon condensation 

B glue « (250MeV) 4 . (5) 

There is an intricate relationship between quark and gluon condensates, as evidenced by 
their behavior as the temperature is increased through chiral restoration T ~ T x sr. Lattice 
gauge calculations || show that about half of the gluon condensate decondenses |5], [II]] as 
the temperature rises through T x sr. Of course, all of the quark condensate^] {qq), which 
can be considered as the order parameter for the chiral restoration transition, goes to zero 
at the critical temperature. 

In the chiral limit (defined as the limit in which bare quark masses are set equal to zero) , 
all of the mass of hadrons such as the nucleon or vector meson can be expressed in terms of 
the gluon condensate Tjf, since the energy-momentum tensor involves only this condensate. 
In the QCD sum rule results, however, the dynamically generated hadron masses depend 
chiefly on the quark condensate {qq), the gluon condensate entering as a correction to this, 
at about the 10% level.0 

The intricate interplay between quark and gluon condensates in the determination of 
the p-meson was made clear by Adami and Brown Jll[ |. The argument given there was 
rather formal, phrased in temperature-dependent QCD sum rules. We present here a much 
simpler physical argument. 

The role of the {qq) quark condensate is to produce most of the p-meson mass. It can 
be considered as generating the mass by producing a repulsive scalar potential, equal in 
magnitude to the mass m p , from the scalar quark condensate in the vacuum, as shown in 
Fig§ 

Clearly the quark condensate couples repulsively to the p meson, the repulsive scalar 
potential representing the />meson mass. However, from the sign of the coefficient of gluon 
condensate, C G i{T 2 , M 2 ), this (gluon) condensate couples attractively. The dimension-4 
operator G 2 where G^ v is the gluon energy-momentum tensor is scaled with M^ 4 , where Mb 
is the Borel mass. With increasing density and decreasing quark condensate, Mb — ► M B , 

#3 We denote the current quark by q and the constituent quark by Q. 

#4 An analogy is that of driving a car. The car moves because the wheels turn (gluon condensate). However, 
the driver determines the velocity at which he is moving from the speedometer (quark condensate). The 
wheels and speedometer are technically connected (QCD). 
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Figure 1: The role of the quark condensate {qq) in generating the p-meson mass can be 
visualized in terms of a scalar field coupling the p meson to the negative quark condensate 
in the vacuum. This simple picture was developed in [O], in order to relate the NJL model 
to this picture of the vacuum condensate producing the p meson mass by this mean-field 
mechanism. 



so decreasing. Although the gluon condensate enters into determining m p only at the 
~ 10 % level at zero density, at finite density, it would seem to become relatively more 
important because of the decreasing (M^)~ 4 effect. The gluon condensate itself is known 
to change rather little with density. At the point where the repulsive contribution from the 
quark condensate has dropped to the value of the attractive contribution from the gluon 
condensate, the p-meson mass m* would go to zero, before {qq)* goes to zero. It seems 
unreasonable that m* should go to zero before chiral restoration. (There is an exception to 
this in connection with Georgi's vector limit which will be discussed in detail below.) 

Adami and Brown [jnj resolved this difficulty by showing that when the perturba- 
tive (black-body) temperature effects are summed and incorporated into the coefficient 
Cq2{T 2 , M 2 ), then the modified coefficient goes to zero as T — > T x $r. The contribution of 
the gluon condensate to the p-meson mass is always small, ~ 10 % of that from the quark 
condensate. Pictorially this can be represented by Figj2|. Thus the gluon condensate enters 
as a "Lamb shift" correction to the p-meson mass. That is, it is the modification by the 
virtual gluon field in the mass dynamically generated from {qq). 

One can thus understand how the dynamically generated m* goes to zero as {qq)* goes 
to zero. However the scaling (jl|) and (0) - "BR scaling" - that results from splitting of \ 
into a x* an d x', has not been clearly explained. A recent paper by Beane and van Kolck 
[13] suggests a remedy to this deficiency. We describe their development in the next section. 
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Figure 2: Contribution of the gluon condensate to the p- meson mass. 



2 A Dilatation-Invariant Low-Energy Lagrangian 

The starting point of Beane-van Kolck treatment is an old "theorem" of Weinberg 
|14j which states that the full content of quantum field theory (here QCD) is given en- 
tirely by general physical principles like unitarity, cluster decomposition, Lorentz invariance 
etc. together with the assumed internal symmetries. QCD is characterized by the global 
SU(2)l x SU(2)r symmetry, which however is not obvious from the spectrum. Weinberg's 
theorem suggests that an "equivalent" field theory exists, which employs constituent quark 
(quasiparticle) fields that realize chiral symmetry nonlinearly and therefore includes explic- 
itly Goldstone boson fields. Brown and Rho jlj] went further, beginning from the Skyrme 
Lagrangian out of which the fermions (nucleons) arose as chiral solitons. This would imply 
a replacement, at a later stage, of the nucleons made up out of constituent quarks in the 
Beane-van Kolck approach by baryon fields. 

The low-energy effective theory that we wish to study must manifest the same confor- 
mal symmetry possessed by QCD. Yang-Mills theory possesses scale invariance which is 
broken by the anomaly while the scale invariance of the fermionic part of QCD is broken 
by the quark masses. Now the variables of the low-energy Lagrangian are the Goldstone 
bosons (pions) and constituent quarks (or nucleons). Low-energy vector and axial vector 
interactions can be viewed as made up out of correlated pions or in a more compact way 



introduced as hidden gauge bosons [15]. The latter allows us to relate the property of the 
vector mesons in hot and/or dense matter to Georgi's vector symmetry [|17]] as we shall 
discuss later. 

The dilatation invariant chiral Lagrangian supplemented by the trace anomaly that 
Beane and van Kolck arrive at is 
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- l -tr(G, v Gn-V(x) + --- (6) 
Here m is the constituent quark mass (which will be defined more precisely later) and 

V, = + (7) 

\ = ^Z%S-Zd^) (8) 

with 

f = U = e in/u with vr = r-vr(x). (9) 

The trace "Tr" is over the flavor group and "tr" over the color group. Transformation prop- 
erties of £ and of the constituent quark doublet ip under chiral and conformal transformation 
are given in J13[| . What concerns us chiefly here is that a sufficient power of the scalar field x 
must be put into each term so that the operator scaling properties, under dilatation, of each 
term in the Lagrangian are those of QCD, i.e, that the Lagrangian is dilatation- invariant, 
except for the potential term which breaks this invariance. This potential subsumes radia- 
tive corrections of high chiral order and hence can be very complicated, the only condition 
being that it gives precisely the trace anomaly of QCD in terms of the scalar field x- The 
precise form of this potential is neither known nor necessary for our discussion. The basic 
assumption here is that the symmetry broken by the anomaly could be represented by that 
of spontaneous breaking and that the Coleman- Weinberg mechanism "chooses" the vacuum 
value of the scalar field in hot and/ or dense medium. 

In the "derivation" of the low-energy effective Lagrangian from a microscopic Lagrangian 
(say, QCD or a fundamental modeling of QCD), scaling is always broken in obtaining the 
current algebra term jf%Tr (d^Ud^U^). Generally, a regularization is made to incorporate 
high-energy degrees of freedom in low-energy effective theories and a momentum cut-off 
A that delineates the high-energy sector from the low-energy one gets transmuted to the 
physical constant f n , the pion decay constant. The scale breaking through the potential 
V(x) renders this transmutation automatic and more or less unique. 

In the next section we shall review how Beane and van Kolck relate the scalar field x to 
an effective scalar a field by introducing a linear basis 

CT + IT ■ 7T = U X- (10) 

It is in the linear basis that low-energy current algebra is reconciled with high energy Regge 
asymptotics and consequently becomes more relevant for high-T and high-density processes 
we are interested in. 
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As shown in Fig.[T] the scalar a field can be thought of as generating a constituent quark 
mass niQ through 

m Q = g a Q Q (a)/ml (11) 

where the subscript Q stands for constituent quarks of light flavors. Compare this to the 
mean-field model depicted in Fig.||, 

(a) = -2g aQQ (qq)/m 2 a , (12) 

the factor 2 coming from the two light-quark flavors. In terms of the constituent quark 
mass rag one can then write down the vacuum energy (see ||, eq.(3.11)) 



B vn r B 



X SR 



, / A d 3 k r- =- 1 f A d 3 k m 2 Q 3A 4 

12 U w^^-zl wjwh; + s? (13) 



where A is a cut-off set to obtain the known empirical value of the quark condensate 

(0\uu\0) = -(240 MeV) 3 . (14) 

Whereas the E vac of (|i~3|) is equivalent, physically, to the condensation energy in the NJL 
model II which involves the quark condensate (0|un|0) that depends on the cut-off A 
quadratically, eq. (|l3|) , however, diverges only logarithmically with A which means that the 
B x sr is reliably determined in this formulation. 

Determining A from the quark condensate (0), one finds 

B x sr « (140 MeV) 4 m 50 MeV//m 3 , (15) 

just the value of the MIT bag constant. We shall see in the next section that this is just 
the value obtained from the Beane-van Kolck theory. 

We see that in the low-energy sector, dilatation invariance is broken by the development 
of the quark condensate. The scale of this breaking is {B x sr) 1 ^ a ■ In QCD, however, scale 
breaking occurs in the gluon sector at QCD loop level. The magnitude of this scale breaking 
is an order of magnitude larger than the B x sr of eq. (|l~5|) ; namely 



£ 9hie ^ 500MeV//m 3 , (16) 

obtained from the QCD trace anomaly. Although B l J^ e ~ 250 MeV is not so much larger 
than B][g R k, 140 MeV, we know that the magnitude of the gluon condensate determines 
the glueball mass 

M GB w 1.5 - 2 GeV. (17) 
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The QCD sum rules tell us [11] that the gluon condensate is relatively unimportant for the 
masses of hadrons made up out of up and down quarks, entering only through the Lamb- 
shift type effect of Fig.^. What this implies is that the scalar x held in the Lagrangian (|6|) 
has mostly to do with the quark condensate 

B xS R oc x 4 (18) 

rather than with the gluon condensate B g i ue . Just how B g \ ue reduces in light-quark systems 
to B x sr is not well understood. In any event, the x* of Brown and Rho M should be 
identified with the mean field of this smoothed field. This will be discussed in more detail in 
the next section. Here we consider the implication of the restoration of the correct scaling 
properties of the operators in the low-energy chiral Lagrangian. 

The mean-field component of x which we denoted x* previously will be a function of 
temperature and density. As in BJ, we define 

K = Ax*. (19) 

We must stress once again that the x field is here already the smoothed field of the quark 
sector that includes no gluon fluctuation. An equivalent way of making the separation 
is to divide the scalar glueball field ^ associated with the trace anomaly to a "smooth" 



component Xs and a "non-smooth" component Xns, X = Xs + Xns-, and "integrate" out the 
non-smooth component Xns from the effective Lagrangian. What appears in the Beane-van 
Kolck theory is then Xs which is relabeled as simply x- It is the mean-field component of 
this field that scales the quark condensate in the BR scaling. Its fluctuating part cannot in 
general describe a single local degree of freedom since it represents multi-pion excitations 
including the continuum but in hot and dense medium can be interpolated by a local scalar 
effective field denoted a. We will have more to say on this point later. 
Similarly we define 

m Q = m QX* (20) 

from which we see that the constituent quark mass vriq scales with temperature and density 
in the same way as the order parameter /*. In particular, m,g must go to zero with chiral 
restoration, /* — ► 0. 

In order to obtain the scaling of the vector meson masses, it is convenient to look as in 
at the quartic Skyrme term 

1 



32 5 - 



r Tr [tfdM,tfd v U] 2 . (21) 



* 5 Instead of dividing the field x, one could choose to divide the trace of the energy-momentum tensor 
~ TrG 2 (where G^v is the gluon field tensor) into a "quarkish" component H q and a gluonic component H g 
as in jL6|. This separation may lead to a different result on scaling, although we have not investigated this 
possibility. 
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One can think of this arising when the vector mesons are integrated out from the effective 
Lagrangian. The constant g can be identified as the hidden-gauge coupling constant. Now 
this term is scale- invariant as it is, therefore no x field need be multiplied to it. When the 
scale invariance is spontaneously broken, the coefficient of this term remains independent 
of the factor x*> s ° we have 

9* = 9- (22) 

As shown recently by Harada, Kugo and Yamawaki flT8| l, the effective chiral theory, when 
hidden gauge symmetry of the vector mesons is implemented, has the exact low-energy 
theorem 

ml = 2g 2 fl (23) 

which is essentially the KSRF relation. This theorem is proven in zero-temperature and 
zero-density regime but we see no reason why it should not hold in medium. This would 
imply 

m*\ = 2g* 2 f*l (24) 

from which, with (|22|), follows 



^ « £ . (25) 
my f w 

The question of how the nucleon mass scales with temperature and density is a lot 
more subtle, because at finite density, a vector mean-field can develop and build up. There 
have been several studies on this issue |22[ but no satisfactory answers have been obtained. 
We postpone this question until after we discuss the effective a field that results once the 
Lagrangian (|6|) is transformed to a linear basis. 

We do not discuss the pion mass, which results from explicit chiral symmetry breaking, 
a phenomenon associated with the electroweak scale which is much higher than the chiral 
scale we are dealing with. From this point of view, the most reasonable thing to do is 
to assume that the pion mass does not scale within the range of temperature and density 
involved in the chiral Lagrangian. 

In 11, we deduced the scaling of (qq) with density from the assumed operator transforma- 
tion properties of the explicit chiral symmetry breaking term in the low-energy Lagrangian. 
The result was 

" <??>*! 1/3 ^ ft 



(26) 

Jtt 

For the same reason as for the pion mass, this scaling relation does not follow immediately 
(without some strong assumptions) from the low-energy effective Lagrangian and in fact 
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may not be quite correct as we shall argue later in connection with the modeling of lattice 
results. 

It should be remarked that the Brown-Rho consideration was at mean-field level. The 
question as to what that corresponds to in the sense of Weinberg's "theorem" has not yet 
been addressed. In a broad sense of effective theories, the BR scaling is to represent the 
tree order in chiral perturbation expansion (more on this later). Higher-order corrections 
will surely modify the scaling relation. For example, from the Goldberger-Treiman relation, 
we find that 

9ttNN = 9A 
m N W 

As we shall review in the next section, the in medium drops from 1.26 to ~ 1 as density 
increases from zero to nuclear matter density pq. The cause of this change has to do with 
the role of short-ranged interactions between baryons that involve multi- Fermi interactions 
absent in the matter-free space. This essentially cancels the change in f n . Thus g^NN /tun is 
expected to remain roughly constant up to nuclear matter density, even though the Brown- 
Rho scaling would have drop, but g^NN remain unchanged. Once g\ has gone to ~ 1, 
however, one would not expect it to change much further, and the BR scaling should take 
over at higher densities. 

Before proceeding, we should address the question as to what the precise role of the 
effective (density-dependent and temperature-dependent) masses and coupling constants is 
in confronting physical observables. As in all field theory problems, those renormalized 
quantities are parameters of the theory that have no physical meaning independently of the 
observables one is discussing. For instance, the m*'s are not by themselves physical observ- 
ables. In some special processes, they could be associated with quasiparticle parameters, 
that is, they can be taken as a pole of the particle Green's function but in general, the pole 
of the Green's function may not be approximated by the m*'s. Residual interactions will 
modify the residues and pole positions. In this respect, it is important to realize that the 
m*'s defined in one theory need not be the same as the m*'s defined in another theory. 
Thus it makes no sense for someone to take a particular hadronic model to calculate what 
he/she defines as m*'s and compare with another person's m*'s calculated in a different 
hadronic model. They can only compare physical amplitudes computed with the model 
wherein the parameters m*'s appear. This also means that a particular scaling relation 
present at mean-field level of one person's theory need not be reproduced in another per- 
son's mean-field theory. As we shall note below, it is a particular reparametrization of the 
fields that gives rise to a simply scaling theory; other reparametrizations need not give the 
same scaling of the parameters although the physics can be identical. In fact, unless one 
defines how to compute corrections, the density (and temperature) dependence of the ef- 
fective parameters with a particular model has no meaning. This point seems to be largely 



9 



overlooked in the field. 

This then raises the question: In what sense the BR scaling is to be understood? 

To answer this question, we recall the principal assumption made in Q - and fur- 
ther elaborated in |2(| - that at each density and/or temperature, we have an effective 
Lagrangian that has the same chiral and conformal symmetry as in free-space with the 
parameters of the theory density /temperature dependent. This can be considered as an 
application of Nambu's "theorem" that whenever spontaneous symmetry breaking is in- 
volved, be that in condensed matter, nuclear or elementary particles, the physics involved 
is a generic a model characterized by different length scales It is not obvious that this 
is always possible, so in our case that is a strong assumption. However we propose one way 
of justifying this approach and it is the notion that the nuclear matter is a chiral Fermi 
liquid as proposed by Lynn [19]. 

It is fairly well-established empirically that the interactions of mesons and baryons inside 
nuclear matter are governed mainly by chiral symmetry. But it is also well-established that 
nuclear matter - and nuclear ground states - cannot be described starting from chirally 
symmetric Lagrangians, at least at low orders. So the ground-state (mean-field) properties 
which do not follow naturally from chiral symmetry and the excitation (fluctuation) prop- 
erties which do are not compatible. These can be reconciled, however, if the nucleus is a 
soliton, that is, a chiral Fermi liquid drop with nuclear matter being a chiral Fermi liquid. 
The chiral Fermi liquid is a soliton solution of the chiral effective theory at quantum level. 
The parameters so fixed in this soliton structure can be identified with the m*'s etc. (In 
the BR scaling, it is the potential V(x) which is supposed to fix the parameters, but we 
do not really know how to calculate it. Lynn's chiral liquid approach suggests how to do 
this.) Fluctuations around the soliton would then have the requisite chiral and conformal 
symmetry of the original Lagrangian, so must retain the generic a model form. (We argue 
below that it is a linear a model.) In looking at the scaling of the a field in terms of 
phenomenological meson theories that we describe in the next section, we are implicitly 
assuming this structure. 

While we find the approach described above quite plausible, we must admit that we 
have not yet found a realistic chiral liquid soliton which could be used as the background 
field to build a theory to "determine" the parameters in a self-consistent way. One can 
however make the following conjecture. To the extent that Walecka's mean field theory de- 
scribes nature fairly accurately, the chiral liquid soliton, when correctly formulated, should 
resemble Walecka's mean field theory. The parameters appearing in the effective Lagrangian 
describing fluctuations around the chiral liquid (mean field) will have BR scaling incorpo- 
rated naturally. Indeed we shall see later that this is what is observed phenomenologically 
when we consider fluctuations in the strangeness flavor, i.e, in kaon-nuclear processes. For- 
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mulating this in a rigorous way remains an open problem. 



3 Transformation to a Linear Basis; Interpretation of the a 
Field 

Beane and van Kolck transform their Lagrangian to a linear Lagrangian by making the 
field redefinition 

£ = U X (28) 

where 

S = CT + ir-7r. (29) 

(Note that since the Sugawara field U = e m ^ w is conformally invariant, the tt must scale 
as f n and consequently we have a new tt* for each temperature or density. We shall not 
make this distinction here.) The Beane-van Kolck Lagrangian obtained by the above field 
redefinition is rather involved, but it simplifies considerably in the so-called "dilaton" limit 
(m CT 0, g A ^ 1) to 

C = iQPQ-^tv(G^G^) + 19^-8^ 

j ffi 

+-d ti ad^a - —Q[cr - i-y 5 r ■ ir}Q 

^ J IT 

+ + - 2 ) 2 - ^ 2 + - 2 ) 2 ln[(^ 2 + vr 2 )/^] + . . . (30) 



Since in the vacuum the pion field has zero VEV and we neglect pion fluctuation, we can 
identify 



"4 4 ^4/ 2 
16/1 16 

= (0|a|0) 



since <ro = f n in the vacuum. We find immediately that 

B%~\^f^U- (32) 

For B^g R ~ 140 MeV, m a ~ 918 MeV, in the range of a masses in the bare vacuum 
consistent with the Weinberg- Tomozawa relation. This B x sr is similar to the condensation 



energy in the NJL model, eq.(15), but has the advantage that it follows in a straightforward 



way from the development. Furthermore, no explicit cut-off is needed in this way of arriving 
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Figure 3: Uncorrelated two-pion exchange involving nucleon and A intermediate states 



at the condensate energy, with the bag constant arising naturally at the tree level. So, 
although we believe the NJL model captures the essential physics of the low-energy scale 
bag constant, the Beane-van Kolck treatment is highly more preferable and appealing. 

Let us now discuss the a field which has emerged in going from the non-linear realization 
of chiral symmetry in the transformation to the linear basis. We should mention to start with 
that it is not a mere field redefinition which would leave physical quantities unchanged. In 
going from the nonlinear structure to the linear structure, a nontrivial physics information 
has been injected. This point was explained by Adami and Brown || using a different 
argument. The crucial point is that in the transformation, one has chosen the component 
of the x field such that the a is now composed mainly of even powers of the pion field tt in 
the expansion of the U = e tn ^^ field. The fluctuating part of this a therefore interpolates 
two-pion excitations with a mixing to 4-7T, 6ir etc. The physical implication of this structure 
can be appreciated by recalling the origin of the scalar attraction in the nucleon-nucleon 
interaction. 

Beginning from the NN — > 2ir helicity amplitude, which can be obtained by analytic 
continuation from 7riV scattering amplitude, the scalar attraction can be built up |2j|. A 
dynamical reconstruction of these results produces the scalar exchange in the Bonn potential 
p4| . In addition to the six uncorrelated two-pion exchange diagrams shown in Fig.||, there 
are the contributions in which the two pions scatter on each other in intermediate states, 
shown in Fig.|I]. 

The tt — tt rescattering necessary to form the correlated two-pion exchange is, of course, 
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Figure 4: Correlated two-pion exchange involving nucleon and A intermediate states 



automatically included in the dispersion theory formalism, where the NN — ► 2ir helicity 
amplitudes are obtained by analytic continuation of the irN — ► irN scattering. However it is 
more convenient for our purpose to use a dynamical model [^] since we shall later introduce 
effects of finite density and/or temperature. The correlated two-pion exchange provides 
about 2/3 of the total 2ir exchange in the exchange of scalar degrees of freedom. The 
mass distribution is broad, but the correlated exchange can be handled rather accurately 
by replacement by a sharp mass particle a'. We shall later drop the prime on this object, 
because we claim that this is just the low-mass dilaton of Beane and van Kolck. Effective 
coupling constants and masses for this effective scalar particle are given in various versions 
of the Bonn potential. 

The main rescattering in the Durso- Jackson- Verwest work [25] comes from crossed- 
channel p-meson exchange. Thus, we can picture the a' as resulting from the summation 
of the diagrams shown in Fig.|5[ 

We are now ready to investigate the density dependence of the a mass. In eq. (|25|) , we 
showed that from the KSRF relation mj, /my ~ i e, the vector meson mass drops as 

the order parameter decreases. This density-dependent />meson mass has been included in 
the integral equation, shown in Fig|| determining the a From the NN — > 2tt helicity 
amplitudes supplied to us by the authors of this work, we can track the downward movement 
in mass of the (distributed) scalar strength as the density increases. Since the downward 
shift of the scalar strength originates from the density dependence of the p-meson mass, it 
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Figure 5: Description of the scalar meson employed in terms of correlated two-pion exchange. 
The two pions interact chiefly through p-meson exchange in the crossed channel. 



is clear that for low densities where the term linear in density suffices, that 



m 



* 



« (33) 
m a my 

The phase of the f+ =0 (t) helicity amplitude goes through tt/2 for \ft « 500 MeV, 
showing that the fictitious scalar particle of mass « 600 MeV used in boson-exchange 
models to mimic the enhancement in the /+ =0 (i) helicity amplitudes for NN scattering 
becomes a real resonance of mass « 500 MeV by p « po. 

Thus, for density p ~ po there is a resonance in the scalar channel at m a rs 500 MeV. 
This light mass cr suggests that chiral symmetry be realized in the linear u-model as a 
manifestation of "mended symmetry" as emphasized by Beane and van Kolck. While in 
the zero-temperature and zero-density regime, Nature prefers the nonlinear realization of 
chiral symmetry for which we have now rather compelling evidence, as density/temperature 
is increased in hadronic matter, the dilaton degree of freedom, frozen in the matter-free 
vacuum, gets identified with the effective scalar channel that is interpolated by a single 
scalar a; at p w po the nonlinear realization of chiral symmetry cedes to a linear realization 
with a low-mass a "seen" in boson-exchange models of the nucleon-nucleon interaction. An 
important point in our proposition is that in medium, chiral symmetry is still manifest in 
Goldstone mode up to the chiral phase transition and the dynamics can still be described 
by a chiral Lagrangian. This is in line with our proposition that Lynn's chiral Fermi 
liquid (which we tend to identify with Walecka's mean field) is a background around which 
fluctuations be made. 

The coupling constants furnish an interesting interplay of the "vacuum structure" dis- 
cussed up to now and short-range nuclear correlations that figure in nuclear many-body 
systems. At zero density, the pion coupling constant at tree level is known to be increased 
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by loop corrections, say, from g^ = 1 to qa = 1-26; i.e, 



g-KNN = {g-KNN)treegA- (34) 

Sometime ago, Rho f p7| l and Ohta and Wakamatsu [^] independently predicted a strong 
density dependence in g& 

9A{p) r 8 ^jva 2 (gp)jVA r i 

9A(0) 1 9 l ) P lor 1 

= [1 + bp/po]- 1 (35) 

where ujr ~ 2.\m- K is the A./V mass difference and 

6 « 0.8(^)tva. (36) 

Here (g'^NA is the Landau-Migdal local field correction in the channel AW <-> ATA which 
can be interpreted as a four-Fermi interaction in the effective Lagrangian that figures in 



many-baryon systems |2C]. The Landau-Migdal parameter is usually taken to be ~ 1/3; 



with this value we expectp 5 



9A(po) ~ 1. (37) 

This rapid drop in g& means that in nuclear matter calculations, it is better to use the 
pion-nucleon coupling given at tree level, 

{gTrNN)tree _ 1 / 9-kNN \2 ^ g g /gg-v 

47T _ 4tt 9a 

Since the a now forms the fourth component of SU(2) x SU(2) ~ 0(4), i.e, the linear 
realization of chiral symmetry, it should have the same coupling constant (1381). Indeed this 



is essentially the coupling of the a in the Bonn potential [31] where the versions A, B and C 
of the relativistic OBEP have g'o- J 47r = 8.8, 8.9 and 8.6, respectively. We thus come to the 
conclusion that at p ~ po, there is a satisfactory linear realization of chiral symmetry, with 
7T and a mesons coupled with approximately the tree-level pion coupling to the nucleon. 
(Note that the in-medium pion-nucleon vertex in the 75 coupling - which is what enters 

* 6 Careful analyses of Gamow- Teller transitions in nuclei reveal that the effective qa in nuclear (5 decay is 
(modulo shell effects) near unity already in light nuclei where density is not high enough for the Rho-Ohta- 
Wakamatsu (ROW) mechanism to be fully effective. This can be understood as follows. In light nuclei, it 
is the core polarization effects induced by tensor forces that play the main role in quenching the <?a, 
with the ROW mechanism being less effective. As density increases in heavy nuclei, the tensor force gets 
suppressed by the BR scaling as shown in [^o| and hence becomes less effective in the quenching while the 
ROW builds up. In some sense, the tensor correlation and the ROW mechanism play a complementary role 
and both seem related - albeit indirectly - to the chiral restoration phenomenon. 
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in the linear a model - is effectively of the form, (gnNN)tree/'m i k ~ 9nNN/rnN, so as long 
as the pion mass does not scale, the pion-exchange interaction remains unmodified. This 
feature of the linearly realized chiral symmetry is consistent with the BR scaling.) This 
indicates that nuclear matter manifests a "mended symmetry" envisaged by Weinberg |32|] 
in the large N c limit of QCD. 

The density dependences of the type (|35| ) are normally not included in the nuclear matter 
calculations, so it must be somewhat coincidental that the coupling constant, essentially 
that of eq. (|38D , can be used at zero density in the Bonn potential to describe the nucleon- 
nucleon scattering. At zero density, the low-mass a is not developed, and the full panoply of 
uncorrelated and correlated two-pion exchange had to be worked out. We note that in going 
from p = pq to p = there are two opposing tendencies. The first is that the resonance in 
the correlated two-pion exchange moves up from ~ 500 MeV at p ~ po to ~ 1 GeV [i.e, 
/o(975)] at p = 0. The second is that the coupling increases from (gl NN /4ir)t re e ~ 8.88 to 



(.9aNN 7^71") ~ 14. Even so, Durso, Kim and Wambach [26| find that the scalar exchange at 
p = po is substantially more attractive than that at p = 0. Brown and Machleidt [33] show 
that density-dependent loop corrections must be introduced in order to achieve saturation 
of nuclear matter. If these higher-order corrections are introduced into the calculation of 
m a(po), h w ih not come out as low as 500 MeV. This work is in progress. 

The Bonn potential is undergoing some modifications at low momentum scales resulting 
from imposition of chiral constraints on the tttt scattering, which force the off-shell amplitude 
to go repulsive |34j]. In an earlier work, without these constraints, Schuck et al |35j found 
that attractive tttt interactions produced very strong effects in the scalar- isoscalar (a- meson) 
channel, correlating an unreasonable amount of S-wave strength in the region close to 



the threshold at 2m n . Aouissat et al [34| show that when the tttt scattering amplitudes 
are constrained by the chiral symmetry a la linear cr-model |3(| the strong near-threshold 
strength disappears. 

In the linear <r-model, from the potential 

V = U(- 2 + - 2 ) 2 -f 2 S (39) 



4 

the a mass is 

ml = 2Xfl (40) 

We do not expect the constant A to change with density or temperature, so ml must scale 
with /*: 

< « (41) 
This is the same scaling as obtained in eq.(p5|) for the vector mesons. 
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Once the low-energy tttt scattering is set by the linear <r-model the introduction 
of the /9-meson between pions, as discussed above, must be accompanied by subtractions, 
so that the S-wave irir scattering phase shifts are not destroyed. Built in in this way, the 
/>meson plays a relatively minor role, unlike that in the older work, in the composition of 
the effective c-meson. 

We should mention that the bag constant B x sr of eq.(|l5|) is small compared with B g i ue , 
but large in comparison with energies in low-energy nuclear physics. In the nuclear many- 
body system, it represents an energy of ~ 50 MeV//m 3 . At nuclear matter density, the 
density of the nucleon is 1/6 per /m 3 , so that the nucleon rest mass energy is ~ 150 
MeV//m 3 , and the net binding energy per nucleon of 16 MeV is only ~ 3 MeV//m 3 . Thus 
B X SR is actually very appreciable. 



4 Chiral Restoration at Finite Temperature 

Lattice gauge calculations are now sufficiently accurate to make statements about finite 
temperature behavior of various quantities, as we discuss. The situation is not so good with 
finite density lattice calculations, so we shall not pursue this matter here. 

We first put the finite-temperature lattice results into physical units following Koch and 
Brown || who used the results of Kogut et al [|| calculated with 6 time slices. Later, more 
extensive results with 8 time slices [[37]] are found to more or less confirm the 6-time-slice 
results. 

In lattice gauge calculations, the quark condensate (qq) is measured for each tempera- 
ture. In an attempt to understand the lattice results in terms of effective fields, the entropy 
was computed in [38] as a function of temperature, assuming the increase in entropy to 



result from the heavier hadrons, p, uo and a±, going massless. The number of degrees of 
freedom in hadrons was, however, limited to 24, the number of quarks, since the latter are 
the fundamental degrees of freedom. In a modeling of the lattice data along the same line, 
Koch and Brown || calculated the rate of increase of entropy in the quarks under the two 
scaling assumptions: 



1. That the hadron masses scale to zero as 

w 1/3 



(42) 



# 7 Wc are aware of the fact that the linear cr-model without matter coupling is not consistent with chiral 
perturbation theory. Here we are focusing on low-energy S-wave -kit scattering in the presence of matter 
fields such as nucleons and vector mesons, so that the linear realization of chiral symmetry does not in 
practice suffer from this fundamental defect. 
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2. That they scale as 



m h (qq) 



(43) 



The scaling (42) is what was obtained in ^ with what we now believe to be too naive an 
assumption ^ while the scaling d43| ) is obtained in the NJL model for the scalar field a. 
The NJL model cannot make predictions for vector mesons whose mass scale is of the same 
order as the cut-off, so the model is moot on the p and a±. 

The scaling (|4^) which we shall call "Nambu" scaling gives an excellent fit to the lattice 
results, the lattice entropy increasing with temperature substantially faster than that cal- 
culated from the scaling (^2|). However, one should be cautious in interpreting the results 
as a support for the Nambu scaling. The argument that goes into this theory, namely that 
the first 40 (24 quark and 16 gluons) degrees of freedom go massless and nothing happens 
to the masses of the other particles, is extremely crude. It resembles in some sense the 
Debye theory of phonons, where there is a sharp cut-off in the phonon spectra at the point 
where the degrees of freedom in it equal the number of underlying degrees of freedom. It 
must be admitted, however, that crude though it may be, our Debye-like theory, with linear 
dependence of the hadron masses on the quark condensate fits the lattice results for the 
entropy surprisingly well. It should be stressed that this fit uses only the linear relation 
between the hadron masses and the condensate. The temperature comes in only in the 
Boltzmann factor when one calculates the entropy. The temperature may be somewhat 
incorrect because Koch and Brown obtained it from the asymptotic scaling relation, i.e, 
from the lowest-order perturbative expression. The color coupling constant g is not small 
enough for this to be quantitatively correct. Thus the temperature scale in Fig|] may not be 
accurate. Using a better temperature scale will not, however, change the fit of the "Nambu" 
curve to the lattice data appreciably. 

A caveat in what we have been discussing here is that it might not be a good approxi- 
mation to have all the hadron masses scale universally and hence it would be premature to 
rule out one option in favor of the other. What may be considered solid in our argument is 
that the lattice results provide an empirical support for the scaling of the masses together 
with the connection to the quark condensate as given by (|43|). 

Another interesting information one can get from the lattice results is the disappearance 
of the dynamically generated masses of constituent quarks (or hadrons) . Lattice calculations 

# 8 The reasoning based on chiral symmetry and trace anomaly that leads to the BR scaling (|l|) does not 
uniquely lead to this relation. As mentioned below eq.(^6|), there is an additional assumption that goes into 
it, namely, a relation between the quark scalar density and the chiral field U which is not precisely given. 
As such this relation is not to be identified on the same footing as eq.(Q) although we loosely call it a BR 
scaling. 
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Figure 6: The Koch-Brown analytical fit || to (^ip), normalized to unity at T = 0, is 
shown by the dropping dashed line. The data from lattice gauge calculations || is shown 
without error bars. The effects of the bare quark masses in ref.|| were removed by Koch 
and Brown, so that their (dashed) line drops to zero. From the lattice results, the entropy 
can be calculated and the result is plotted, again without error bars, as the rising dashed 
line. The lower solid line follows from the scaling (^); the upper solid line follows the 
"Nambu" scaling of ©. 
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measure screening masses of hadrons. DeTar and Kogut |39| found in four-flavor calculations 
that the screening masses of the p and a\ mesons come together, as expected, at T = T x sr. 
Furthermore, to the accuracy measured, the common screening mass came out to be 

uj scr « 2vrT; (44) 



i.e, irT per quark. Later two-flavor calculations [37] confirmed these results and found 
further that the common screening mass for the nucleon and its chiral (parity-doublet) 
partner was 

fa 3vrT (45) 

for T > T x sr. As shown by Adami and Brown |J, if a residual mass mo remained above 
TySRj the screening mass for the vector mesons would be 

u scr « yl^ + ml « ttT + 1^. (46) 

In fact to the accuracy calculated in lattice gauge simulations, no mo is needed, the conclu- 
sion being that - to this accuracy - mo is zero. Of course a small mo could easily escape 
detection. 

Direct calculation of the constituent quark mass (not the screening mass) for four flavors 
yields a mass which drops to the perturbative value 



m, 



eff=9 



T/V6 (47) 



for T going from T x sr to 1.75T x sr. Similar calculations of the p- meson mass [41] find m* 
to decrease rapidly as T increases. There is of course a caveat to this. One has to be careful 
in applying four-flavor results to the two-flavor world, since the chiral restoration transition 
is first-order in the four-flavor case and the quark condensate does not seem to change much 
below T x sr, whereas the transition is a smooth second-order one in the two-flavor case. 



Calculations of the Bethe-Salpeter wave functions of ir- and /j-mesons [42] show these 
mesons to be more compact, i.e, more closely correlated, at 1.5T x sr than at zero tem- 
perature. Koch et al j£| have shown that this can be simply understood in terms of a 
dimensional reduction of QCD at high temperature. 

The results of ref.jljj could be summarized as follows. The quarks propagating space-like 
experience a space-like potential that rises linearly at large distances, as first calculated on 
the lattice by Manousakis and Polonyi [44]. In a "funny space" obtained by interchanging z 



and t, the quark-antiquark wave function is calculated from the non-relativistic Schrodinger 
equation in which the effective quark mass is 

m eff = TrT. (48) 
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As with the screening masses, were a dynamically generated mass still present above T x sr, 
the effective mass to be used for the quark would be J tt 2 T 2 + mjj where as before niQ is 
the dynamically generated mass. Similar calculations of the it- and p-meson wave function 
were subsequently made by Schramm and Chu |49] who extended them to a wider range 
of temperature. They found that the wave functions changed even less with temperature 
than in the calculation of Koch et al [[f^] . 

There are two remarkable features in the theoretical results of Koch et al: (1) Dimen- 
sional reduction of QCD seems to work well right down to T x $r ~ 140 MeV, a surprisingly 
low energy for something which is supposed to happen at "asymptotic" energy^; (2) the 
helicity-zero state of the p meson comes out degenerate with the pion while the helicity 
± states of the vector meson are degenerate with each other. This is easily understood 
with the quarks massless in the dimensional reduction. In the "funny space" in which z 
replaces T, at least asymptotically in temperature, the configuration space with the new z 
becomes two-dimensional with only x and y directions. Spin must be either perpendicular 
to the (x, y) plane or lie in the plane. The p meson has gone massless and behaves like 
a (charged) photon with helicities ±1 perpendicular to the plane. The helicity-zero state 
which originated from the longitudinally (in plane) polarized component of the p before it 
went massless now behaves as a scalar and forms a multiplet with the charged pion. We pro- 



pose that this situation corresponds to the Georgi vector limit 17] in which the longitudinal 
components of the massive p meson decouple and become degenerate with the pseudoscalar 
Goldstone bosons, the pions, thereby the p becoming massless. Georgi envisaged this limit 
as arising in some particular limit (such as large N c limit) of QCD. It appears from the 
results described above that the Georgi vector limit is relevant in the vicinity of the chiral 
phase transition. Note that this is also in line with (though not equivalent to) Weinberg's 
"mended symmetry" scenario exploited by Beane and van Kolck. p^l We shall have more 
to say on this matter in connection with quark number susceptibility in the next section. 



* Some insight into why the effective dimensional reduction works so well has been given by detailed 
calculations of Suzhou Huang and Marcello Lissia [Q. As noted above, the lowest Matsubara frequency 
becomes the chiral mass in our dimensionally reduced space. Corrections to this dimensional reduction 
involve only the higher Matsubara frequencies, the next one being 2nT higher in energy. They are prefixed 
by the running coupling constant a s . Because of the 2tyT, one might expect the scale A 2 to be used in 
obtaining this a s = a s (A 2 ) to be a factor of ln(27r) 2 higher than that for the a B in the dimensionally reduced 
space. The detailed calculation of Huang show it to be a factor of ~ 80 larger (for N c = Nf =3), so that 
a s is cut down by (ln80) _1 ^ \. 

* 10 S. Beane in private communication pointed out that the Georgi vector limit and Weinberg's mended 
symmetry are basically different, particularly near a phase transition. 
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5 Quark Number Susceptibility; Description of the Vector 
Interaction for T > T x sr 



The quark number susceptibility [47, 48 1 is defined as 



X± = (d/dp u ± 8/dfJLd) ( Pu ± Pd ) (49) 

where the + and — signs define the singlet (isospin zero) and triplet (isospin one) suscep- 
tibilities, jJL u and pa are the chemical potentials of the up and down quarks and 



Pi = Tr Niexp 



i=u,d 



/V=m))lV (50) 



with Ni the quark number operator for flavor i = u.d. Generalizing McLerran's expression 



[49] for the susceptibility and changing his baryon number density to quark density, we 



obtain 

X± = {VkT)- 1 I d A x{{{N u (x) ± N d (x))(N u (0) ± N d (0))}}. (51) 



The x± are called the singlet and triplet susceptibilities. We shall see that the x+ ls m the 
cu-meson channel and the X- m the p-meson channel. We can thus call x± the isosinglet and 
isovector susceptibilities. The susceptibilities measured on lattice are given in Fig.^. For 
<577(2) symmetry, we expect x+ = X- an d this is what one observes in the lattice results. 

Although the singlet susceptibility has large errors, it is statistically consistent with the 
more accurate nonsinglet susceptibility. We shall discuss the latter here. It can be seen 
that the susceptibility is small at low temperatures, rising rapidly as T moves through the 
critical temperature T x sr up to ~ 70 % of the free-quark value, designated x^ i n Fig.^. 

We can distinguish three separate regimes in temperature, which we discuss one by one. 
In the very low temperature regime - in which we are not particularly interested - the 
nonsinglet susceptibility is saturated by the p meson [pq] 



X NS |r«T xSfl , « (/cT) J d A xG pp (x) (52) 

where G pp is the p-meson propagator. 

As the temperature T moves upwards to the onset of the phase transition, the constituent 
quark model would be more appropriate |50| |. In RPA approximation of this model as 
depicted in Fig.^, the susceptibility below the critical temperature is 

X = Xo/(l +9vXo) (53) 

where g v is the coupling of the constituent quark to the vector meson and xo is the suscep- 
tibility for non-interacting quarks which at T T x sr where the dynamical quark mass rriQ 
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Figure 7: Quark number susceptibilities as calculated by Gottlieb et al p7[ . The singlet 
(isoscalar) susceptibility is given on the left-hand side and, the non-singlet (isovector) on 
the right-hand side. The free-quark susceptibilities Xs\s^ corrected for effects of the finite 
lattice size, are also shown. Horizontal arrows label the values of x m the continuum (where 
X = 0.125 in lattice units for 4 time slices) and corrected for the finite size of the lattice 
(labeled 8 3 x 4). The stars on the figure show our values calculated from perturbative gluon 
exchange, eq.(|67|) and equations that follow. The star at T c has to be moved somewhat to 
the right, because the phase transition is a smooth one, of width ~ 20 MeV @ and our 
calculation applies only when the transition is completed. 
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Figure 8: The quark number susceptibility below T x sr is described in RPA approximation 
by summing quark-antiquark bubbles interacting by exchange of p mesons. 



has dropped to zero has the value 

Xo « N f T 2 (54) 

with Nf the number of flavors. Now expressing the constant g v in terms of the vector gauge 
coupling g that figures in the hidden gauge symmetry Lagrangian | 



1 g 

9v ~ ~. — o-< (55) 
Am 2 



v 



Thus 



X(T) ~ • (56) 

l + i£xd(T) 



4 m- 

Kunihiro [|50| investigated in the NJL model what happens to the susceptibility as 
expressed in eq.([56|) as temperature approached T x sr from below and concluded that the 
vector field should decouple to explain the rapid enhancement of the susceptibility observed 
in the lattice results. This means that g 2 /m v must steeply go to zero. The NJL model 
cannot explain this decoupling since the constant g v in that model is a relic of degrees of 
freedom integrated out in arriving at the NJL form of effective theory and hence there is no 
way to know how this constant runs as a function of temperature or density. On the other 
hand the hidden gauge symmetry theory tells us at least qualitatively how the constant 



might run. The hint comes from the recent result of Harada and Yamawaki [51|. 



Using the hidden gauge symmetric Lagrangian [15], Harada and Yamawaki showed with 
one-loop j3 functions that the gauge coupling g = 2(1 + K)g pwn (where k is a parameter 
#n We are not putting asterisks but the temperature dependence of the constants and masses is understood. 



24 



which takes the value k = — 1/2 in free space) scales to zero as g(fi) ~ (hi/x) 1 with k — > 
when fjL — > oo. As the gauge coupling goes to zero, the vector-meson mass which is given by 

m 2 v = (l + K)- l flg 2 (57) 

goes to zero. We should not take this perturbative argument too literally for a quantitative 
understanding since there could very well be some important non-perturbative effects in 
this regime that could modify the running of the coupling constant but our assumption is 
that it is very possible that the gauge coupling constant drops - as suggested by the lattice 
data - faster than logarithmically When g = 0, then the vector mesons decouple, the 
longitudinal component of the vectors becoming scalar Goldstone bosons degenerate with 
the pseudoscalar Goldstone bosons and the vector meson becomes massless. This is the 
Georgi vector limit [17]. In this limit, we will be left with 



X ~ X0. (58) 



Before the vector decoupling leading to (|58| ) , we can use the KSRF relation at T near T x sr 
(which seems to be justified by the work of Harada, Kugo and Yamawaki |l^]) and xo ~ 2T 2 
for T ~ T x sr to get the ratio just before the critical point 



X(T< SR )/ X) (T< SR ) « -j—S - 0.47 (59) 



for T x sr ~ 140 MeV and k = 0. Here we are assuming that as suggested by lattice results, 
/tt remains at its zero temperature value up to near T x sr. (The constant f n is believed to 
fall very rapidly to zero within a small range of AT near the critical temperature.) The 
ratio ( |59| ) is in agreement with the lattice data at T < T x sr. 



#12 The phase with g / preceding the Georgi vector limit contains the scalar Goldstone bosons (s) that are 
the longitudinal components of the massive vector mesons. Since k — with f s = f v , the symmetry SU(2) x 
SU(2) is restored for the would-be scalar Goldstone bosons s and the pions. An interesting observation to 
make here is that when k = 0, the direct photon coupling to the charged pion which is proportional to 
(1 — 2(i+«) ) e becomes | whereas in the normal phase where k — —1/2 the direct coupling vanishes. This 
means that as temperature and/or density is raised, the photon coupling deviates from the canonical vector 
dominance picture. It also means that the photon couples half-and-half directly with the pion and through 
the massive p for j / 0, c = and with the massless quark-antiquark pair of the p quantum number for 
ft = g — 0. It should however be pointed out that as discussed by Pisarski|H[, it is possible to preserve 
the vector dominance at all temperatures. This option would then violate hidden local symmetry with the 
consequence that the p mass will go up - instead of down as in BR scaling - as temperature increases. 
While we favor the hidden gauge symmetry prediction which is consistent with the notion that the vector 
meson mass is dynamically generated, nothing rules out the vector dominance option and it will be up to 
experiments to decide which is chosen by nature. This makes the planned dilepton measurements of the 
mass shift of the vectors at GSI, CEBAF and other laboratories particularly tantalizing. 
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Figure 9: Perturbative calculation of x- The quark and antiquark coupled to the quark 
density p q interact by exchanging gluons depicted by wavy lines. 



Let us finally turn to the third regime, namely above T x $r. It has been shown by 
Prakash and Zahed [52] that with increasing temperature, the susceptibility goes to its 
perturbative value which can be calculated with the perturbative gluon-exchange diagrams 
of Fig.||. The argument is made with the dimensional reduction at asymptotic temperatures, 
but as mentioned above, it seems to apply even at a temperature slightly above T x sr. 

In the following, we schematize the Prakash-Zahed argument. Let us assume, in accor- 
dance with the Georgi vector limit, that the vector meson exchanges are decoupled in this 
regime. Above T c , the gluon exchanges are just those of Koch et al JET f^ 16 - 



with 



V(n) = ^a zA a z , 2 5(r t ) (60) 

e 2 - ~/T (61) 

for QCD (with g the color gauge coupling) and 5(rt) is the (^-function in the two-dimensional 
reduced space. Here m = ttT is the chiral mass of quark or antiquark as explained in [43]. 
Possible constant terms that can contribute to eq.(pO|) will be ignored as in [f43|]. 

In order to evaluate the expectation value of the S(rt), we note that the helicity-zero 
/>meson wave function in two dimensions is well approximated by 

i> p « Ne- rt ' a (62) 

with a ~ | fm and the normalization 

N 2 = -L (63) 
na z 

#13 Note that <r z ,icr Zi 2 = —1 for both the pion and the helicity-zero component of the p which is degenerate 
with it. This degeneracy of the helicity-zero states should be checked directly by lattice gauge calculations. 
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For the helicity ±1 p-mesons, <t z ,i&z,2 = 1 ; so we find that the expectation value of V is 

K > 3<K 2 T 2 a 2 { ' 

Going from the lowest order process xo to the one involving the first rung in the ladder, 
Fig|9|, means introducing the correction factor 

~ 27rT' (65) 

the denominator being the unperturbed energy of the quark-antiquark pair. Summing the 
rungs, as shown in Fig.||, to all orders gives us 

2L = (i + ^T) . (66) 



Xo V 2vrT 

The lattice calculations jl7| use 6/g 2 = 5.32 which would give a s = 0.07 at scale of a -1 
where a is the lattice spacing. (The relevant scale may be more like 2n/a.) Calculations 
use 4 time slices, so the renormalized g is that appropriate to a _1//4 . Very roughly we take 
this into account by multiplying the above a s by ln4 2 ; therefore using a s = 0.19. With this 
a s and the above wave function, we find 

-^ir « 0.68. (67) 
Xo(T c ) 

As can be seen by the * at T c on the right-hand graph of Fig.^, this is just about the ratio 
obtained. 

In view of the crudeness in our determination of a s , this quantitative agreement with 
lattice results may not be taken seriously. However, it should be noted that the perturbative 
correction (V)/27rT goes approximately as T~ 2 , neglecting the (logarithmic) change in a s , 
and it is seen from Fig.^ that this T-dependence fits that of the calculated x/xo quite 
accurately for T > T x sr. 

Going towards the vector limit, chiral SU (2) x SU (2) symmetry is "mended," so 

U = fs (68) 

where /„■ and f s are the constants defined by 

(0|4|^'(g)) = iUq^, (0\V^(q)) = if s q^ (69) 

where V* (AV) is the vector (axial-vector) current. The isovector scalars 5* correspond to 
the longitudinal components of the p. (We are ignoring here the dilaton a discussed above.) 
Thus while for g / 0, ^ ins, the equality (|68|) still holds by the mended symmetry. 
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Now assuming that the matter above T x $r corresponds to the vector limit with g = 0, is 
the relation (|68[) expected to hold ? 

As discussed above, in the dimensionally reduced "funny space," the tt and the helicity- 
zero p are degenerate (with their dynamical mass equal to zero) and their wave functions 



become identical Ji3[. Therefore we do expect (68) to hold trivially. This is of course 
consistent with the picture of the ir and p made up of two non-interacting massless quarks 
with the Matsubara frequency ttT per each quark. A short calculation of the matrix element 



with the 7r wave function in the "funny space" gives, for large T^ 14 



U ~ <VIT (70) 

where c is a constant << 1 and g is the color gauge coupling constant. Of course just as 
the screening mass has no direct physical meaning - though perhaps related to physical 
quantities through analytic continuation, one cannot give a physical meaning to @ in the 
sense of the pion decay constant. This may account for the fact that f n grows with T just 
as the screening mass does, while one expects the physical f n to go to zero as one would for 
the dynamically generated mass denoted above as rriQ. We clearly need to know how to go 
from "funny space" quantities to physical space ones. The link is lacking at the moment. 
What we can say however is that the results of the lattice calculations are consistent with 
Georgi's vector limit encoded in (|6£ 



6 Effects in Heavy-Ion Collisions 

Consequences of the vanishing of the hadronic coupling gy of the hidden local symmetry, 
leaving on the colored gluon exchange, are strong for the relativistic heavy-ion experiments. 
Firstly, it should be noted that freeze-out in the Brookhaven AGS experiments has been 
determined to be [53|^^ 



T f0 = 120 ~ 140 MeV. (71) 

By freeze-out we mean the effective decoupling, in the sense of energy exchange, of pi- 
ons and nucleons. (Less strongly interacting particles, such as the kaons, freeze out at a 

#14 We denote the constant by f v to distinguish it from the physical pion decay constant f n . 

#15 The original determination of T J> 150 MeV from the ratio of isobars to nucleons by Brown, Stachel and 
Welke H was corrected about 10 MeV downward by taking effects such as the finite width of the isobar 
into account. It should be also mentioned that as recently shown in a dilated chiral quark model |5f|, the 
scaling with temperature will be more rapid in dense matter than in matter-free space. Consequently, the 
freeze-out temperature in the AGS experiments where the reaction takes place at a high density must be 
lower than the critical temperature determined by lattice calculations that pertain to zero-density matter. 
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higher temperature, say, T > T x sr.) Lattice gauge calculations give chiral restoration at a 
temperature [56] 



T xSR ^U0 MeV. (72) 

This suggests that freeze-out for particles other than the pion and nucleon is at a tempera- 
ture higher than T x sr and that the pion and nucleon freeze out at about T x $r- This means 
that interactions in the interior of the fireball will be at temperatures greater than T x sr. 
We have already seen that the hadronic vector interaction essentially decouples, leaving 



only perturbative gluon exchange. Hatsuda and Kunihiro [57] show that the pionic and 
scalar degrees of freedom move smoothly through the phase transition and this has been 
verified by the behavior of the relevant screening masses as T passes through T x sr. This is 



not surprising since as shown by Wilczek [58], the linear a model is the Ginzburg-Landau 



effective Lagrangian for the chiral restoration phase transition. It is in this sense that Beane 
and van Kolck |l3| recover the linear a model in the "mended symmetry" regime (and also 
that the BR scaling arguments make sense at tree level). The pion and a fields are just the 
fields of this class of effective theories. 

As noted above, the behavior of the pion mass m n at zero chemical potential (fj, = 0) is 
complicated [[58]]. A common bare quark mass m u ~ = m plays the role of an external 
magnetic field. Basically one falls back onto the Gell-Mann-Oakes-Renner relation 

"4 - CO 

Jn 

extended to finite temperatures and densities. Since both (qq) and f% go to zero (in the 
chiral limit) as T — > T x sr, the behavior of the pion mass is a subtle matter. Indeed for 
T = and low densities, one can use the fact that empirically m n changes only very little 
(if any) with density, which we know from the very small (and repulsive) scalar potential 
from the nucleons in pionic atoms, so we may turn matters around and say that f n scales 
as 

£ocK*7>*| 1/2 (74) 
which would keep the pion mass unsealed [59]. 



The AGS and CERN relativistic heavy ion collisions offer exciting new perspectives 
on this problem. They make it possible to construct an environment of high densities, 
several times nuclear matter density, and high temperature, T > T x sr. As remarked 
above, essentially all of times of the fireball existence, which we shall shortly show to be 
At ~ (25 — 30) fm/c, the temperature is greater than that for chiral restoration. 
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Figure 10: The chiral circle. 



6.1 The Ti-Term Attraction 

Some years ago, Kaplan and Nelson showed that explicit chiral symmetry breaking - 
that gives masses to the Goldstone bosons - could be "rotated away" by condensing the 
Goldstone bosons in dense nuclear medium. Although they discussed kaon condensation, 
which we shall consider later, it is easy to grasp their key idea by looking at the pion. In 
fact, the effect we wish to study is interesting for both the pion and the kaon. Consider 
the "chiral circle," Fig. JO]. At finite density, the energy density from the explicit chiral 
symmetry breaking is given by 

Hx x sb = E* N NN cos 9 + \flml sin 2 6. (75) 

The first term on the right-hand side represents the sum of effects of the explicit chiral 
symmetry breaking contribution to each nucleon mass S^at and the second term the sum of 
the pion mass, with the pion field being given by ir = f n s'm9. Note that the physical pion 
results from small fluctuation about 6 = on the chiral circle. Expanding Hx x SR i n 

H XX SR = const. + -m n I 1 - J fj +■■■ (76) 

where the ellipsis stands for higher orders in 9. We see that the pion has developed an 
effective mas sp lt] | 



* 16 As mentioned above and to be explained below, this effective mass is not "effective" at zero temperature 
because of energy-dependent repulsive interactions that enter at the same order of chiral perturbation theory 
canceling the sigma-term attraction but becomes observable asT-> T x sr. 
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Figure 11: Pauli blocking. 



It is, however, known from pionic atoms that the scalar interaction felt by the pion is very 
slightly repulsive in nuclei. This repulsion is brought about fnl] by Pauli blocking in the 



virtual pair terms, of the type shown in Fig.ll. In pion-nucleon scattering, a repulsion 
develops because of Pauli blocking when the nucleon in the NN bubble in the pion self- 



energy tries to go into the state already occupied by the original nucleon (see Fig. 11 a). With 
use of the Wick theorem, this is described by Fig. lib as a Z-diagram, with the backward- 
going line representing the antiparticle. 

To fully take into account the effect of this type, one has to include also the anti-A for 
pion-nuclear scattering and the anti-decuplet for kaon-nuclear scattering. The virtual pair 



terms of Fig.ll and their decuplet counterparts effectively cancel the attraction expressed 
in the dropping m*, eq. (|77| ) J^ 1 ^ 7 ] 

Up to nuclear matter densities, the A-nucleon energy difference does not change appre- 
ciably. In the Appendix of ref.|^6[, this is explained in terms of the local field correction 
splitting the nucleon and A at finite density. This correction, however, depends upon the 
coupling of vector mesons, p and u, and as the vector mesons decouple at temperature 
T > T x sr, the local field correction should go to zero. Therefore taking the decuplet to be 
degenerate with the octet should be literally correct at high temperatures. 

Our above discussion can be carried over to the case of kaons, where the effects are 



* 17 In the chiral Lagrangian approach of refs,[[5l|, p2[ |(S3| , p4j |, this repulsive contribution is a part of the 
0(Q 2 ) terms proportional to the square of the meson frequency that appear as "counter terms." They are of 
the same order in the chiral counting as the E term of Kaplan and Nelson fed] . This contribution is important 
for both pion-nuclear and kaon-nuclear scattering but as shown in |64[ , it plays an insignificant role for 
kaon condensation in particular. The point of the discussion given here and in the following is that as argued 
in |53|, |(34j this term can be saturated by what we call "pair terms" involving the octet and decuplet whereas 
in general or more specifically in chiral Lagrangian approaches such a term can arise from complicated 
(uncalculable) sources. 
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larger and the situation is a lot more interesting. 
The K + meson is 

\K + ) = \us) (78) 

containing a nonstrange quark and a strange antiquark. The vector interaction between a 
i^ + -meson and the nucleon is, therefore, repulsive and in dense matter 

V k +n = \v NN 90 MeV s -V K - N (79) 
3 Po 

where po is nuclear matter density and Vnn is the vector mean- field potential felt by 
a nucleon. The Kaplan-Nelson scalar attraction, comparable to that giving the pion its 
effective mass, eq.(fT7|), is^^ 



S ™<*f>g-64MeV* 
2mKJ Po 



ik + n ~ • , 2 = "64 MeV ?L = S K - N (80) 



where p s is the scalar density. See also Appendix. At zero density, virtual pair corrections, 



of the type discussed for pions, remove about 37% of this attraction [61]. 

We note that the chief dependence in the S K + N of eq.(|80|) is expected to come, at least 
initially, from that of /. Up to nuclear matter density the decrease in /* is given by eq.(f74|) 
with 

[qq) f*rr4 
Now if we take X^a? ~ 46 MeV, then 

ff ~ 0.6/2 (82) 

thereby increasing the scalar attraction at p = po by a factor of ~ 1.6. Whereas the Sj^+n 
of eq. (p0|) is decreased ~ 37% by virtual pair correction 64], this factor of ~ 1.6 increase 



makes up for them, so we believe that the Kaplan-Nelson term ( |80| ) to be our best estimate 
for nuclear matter density pq. Extrapolation to higher densities is then made linearly from 
here. 

It should be noted that the (repulsive) virtual pair correction factor goes as F = (1 — 
.2>7uj\/m 2 K ). In the case of the K~-meson, the u 2 K decreases with density, the correction 
dropping out. Taking into account the decrease in (/*) 2 of (|82|) , we then find that S^-n 
becomes roughly equal to Vr-n- 



#18 We are using Ti K n ~ 2.83m^- obtained at tree level in pS, but see Appendix. 
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6.2 Subthreshold Kaon Production 



Recent experiments on subthreshold kaon production in nucleus-nucleus collision at ~ 1 
GeV/nucleon in SIS at GSI |)7j give a factor of ~ 3 more K + -mesons than predicted by 
conventional theories |5^] with the best nuclear matter EOS - one with a conventional 
compression modulus and momentum dependence, as we discuss in more detail below. 

Implementing ( |30|) as a correction to the kaon mass, with S^tv = 2.5m,r, slightly less 
than the T,kn = 2.83?^ obtained in |64j], Fang et al |69| manage to explain the enhanced 
kaon production. Similar results have been obtained by Maruyama et al JF9]. In this second 
paper, it appears that the GSI subthreshold K + events are reproduced without decreasing 
the kaon mass with density. But the authors employ the same scalar mean field for the 
A-particle as for the nucleon. However the coupling to the A should be only 2/3 of that 
to the nucleon, the scalar field not coupling to the strange quark. In terms of lowering 
the in-medium threshold* 19 for kaon production, employing the same scalar field for the A 



as for the nucleon is, at least roughly, equivalent to using 2/3 of it for the A and 1/3 for 
the kaon. Now 1/3 of the scalar field, for nuclear matter density p = po, is > 100 MeV 
in magnitude, larger than the Kaplan-Nelson attraction (|80|) & M . So the physics is roughly 
equivalent. 

Fang et al |3^] do not include the (repulsive) virtual pair correction, nor the scaling 
in (f*) 2 of eq.(p^). Introducing these for p ~ po, they roughly cancel. For K + mesons, 



u>k does not change much with density. Initially [64], it increases slightly. Thus, the factor 
F = (1 — 0.37u^/m^) always cuts the S K + N down somewhat. On the other hand, the 
(/*) 2 continues to decrease as p exceeds po, so that m K + will come back down to rriK at 
p ~ (2 — 3) po, the densities relevant for subthreshold K + production. Thus, the scalar 
attraction and vector repulsion roughly cancel at these densities. This seems to be what is 
required to reproduce the subthreshold JC + -mesons[]69|, [7(|. 

We remind the reader that chiral Lagrangians do not admit a scalar exchange between 
the kaon and the nucleon: Scalar fields have no role in chiral expansion. However, in 
phenomenological fits to K- nucleon scattering [72], which do not include the Kaplan-Nelson 



term, an attraction of about the same magnitude [61] as would come from this term is 



introduced by an effective scalar exchange with a coupling constant combination 

9oN N 9oKK / ~ 0.9. (83) 

The K- nucleon scattering is, however, not so sensitive to this attraction as is the sub- 
threshold kaon production, where an attractive scalar interaction substantially increases 

# 19 Randrup and Ko jnj showed that the subthreshold kaon production was chiefly determined by the 
maximum kaon momentum P ma x which in turn is determined by the threshold, given the input energy. 

# 20 With Brown- Rho scaling and our Ekjv, the Kaplan-Nelson attraction will grow up to 80 MeV for 
P = Po. 
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Figure 12: Production of K + and A in NN and AiV collisions from pion exchange. 



the production of kaons. 

From the above discussion, we see that it should be possible to explain enough sub- 
threshold K + production by using attractive scalar mean fields for each u or d quark, 
also the u-quark in the .ff + -meson. This is essentially what Miiller-Groeling et al did, 
although we see from eq.(|83|) that the g a NN9aKK they used is somewhat less than the 
3 (.ffaNN 1 4f) ~ 1-5-2 usually employed in nucleon mean fields. Our discussion of Maruyama 
et al |6^] shows that at least in the determination of P ma x of the K + , this is essentially 
what they did. 

Note that the vector mean field should be coupled to each u or d quark, since it couples 
to the baryon number; thus each u or d quark experiences a vector mean field equal to 
^Vnn(p), as noted earlier. To the extent that the Kaplan-Nelson term (|80l) is equivalent 
to the scalar mean field of ^Snn(p) at quark level, we have the simple picture of the mean 
fields familiar from, say, Walecka theory, being applied at (constituent) quark level. The 
enhanced subthreshold K + production then follows from the scalar mean field essentially 
canceling the repulsive vector one. The vector mean fields do not affect the determination 
of Pmax of the K + , since the baryon number they couple to is the same before and after 
interaction. 

Miskoviec et al |67| make no comparison of their results with the theoretical results 
mentioned above, because of uncertainties in the nuclear matter EOS, and in the AN — > 
KAN and AA — ► KAN cross sections, which produce ~ 85% of the subthreshold kaons. 



Questions concerning the EOS have been sorted out in the past few years |73] and the 
conclusion is that a conventional compression modulus Kq = 210 ± 30 MeV together 
with momentum dependence (as follows from a nucleon effective mass), should be used. 
The elementary processes AiV — * KAN and AA — * KAN have not been measured 
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experimentally. However, the NN — > KAN cross section is given mainly by the process of 
Figjl^. The A's that are present act as a reservoir of energy and the pion in the process of 
Fig, 12 b is more nearly on-shell than in Fig. 12 a. The crucial ratio of AttN to NttN coupling 
is known to be ~ 2 from analyses of many experiments. Thus scaling in this way and 
averaging over the various charge states of the A, as done by Randrup and Ko Iff!]] , should 
be sufficiently accurate to show that there is a real discrepancy between the experimental 
results Ip^l and the theoretical ones [ 



6g [70] . 



The A A?" — > NKA and A A — ► NKA cross sections are being recalculated in Jiilich [75]. 
The strong interaction calculations, such as those of Ko and collaborators mentioned 



above, have been criticized recently by Schaffner et al (76] and Maruyama et al |77|| . Firstly, 
it is pointed out that in calculations such as those of Fang et al [69], the scalar density 
p s = (NN) was "incorrectly" replaced by the vector density py. It is clear that the scalar 
density decreases relative to the vector density as m* N decreases, going to zero as — ► 0. 
On the other hand, the ratio (NN)/f 2 appears in (jS0|), and if the vector mass is brought 
down in medium by the Kaplan-Nelson term, then so is / since, as noted earlier, they are 
related by the KSRF relation m 2 = 2f 2 g v and gy does not change at tree level. Thus, the 
(f*)~ 2 would be expected to more than counterbalance the decrease in (NN) as compared 
with py. This presumption is supported by ref.[7C] which mistakenly applied the same 
scalar mean field to the A as to the nucleon. As noted, for subthreshold K + production, 
this amounts to using 1/3 of the scalar field as the nucleon for the K + ; furthermore, this 
is a somewhat stronger attraction than would be given by the Kaplan-Nelson term without 
scaling the f~ 2 . 

Schaffner et al 76] further suggest associating the square of the vector potential with 
an increase in the kaon effective mass. This quadratic vector interaction was discarded by 
Brown et al [78] on the ground that it meant using, to quadratic order, a term which had 
been derived in the chiral expansion only to linear order. Furthermore these authors argued 
that effective masses should involve scalar fields and that the vector mean field should give 
only a shift in the chemical potential, as in Walecka mean field theory. 

These problems have been simply resolved in the calculation of kaon condensation 

M| where heavy-fermion chiral perturbation theory 



in dense stellar matter [61, 62, |6c 
(HFChPT) was used. In HFChPT, the expansion is made with the velocity-dependent 
positive-energy baryon field B v = e imBl " v ^(l + 7 • v)B. In 64], all terms up to 0(Q 3 ) in 
the chiral counting are summed. The square of the vector potential does not enter at this 
order ^ 2T . As explained in a note added in [64], the difference 5 = p s — py appears at 0(Q 4 ) 
in the chiral counting and taking into account that difference (and also the term of the 

# 21 Formally such a term corresponds to a two-derivative four-Fermi interaction, so that it would come in 
at 0(Q 4 ), one order higher than the order calculated in m]. 
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form of the square of the vector mean field) in chiral perturbation theory as developed in 
would require, for consistency, including 0(Q 4 ) counter terms and calculating leading 



two-loop diagrams. A partial account of such terms as is done in [76, 77 1 would violate 
chiral Ward identities. We claim that the result of |64j| corroborates this point: The large 
effects claimed by Schaffner et al |58| and Maruyama et al [69] do not appear, in that the 
scalar field here plays a relatively minor role. Its magnitude is varied more than a factor 
of two between calculations in and [64], but the results are not so different. It thus 
appears that the schematization of [68] and |3(| to relativistic mean field theory is too crude 
an approximation. We also point out that the argument of Schaffner et al |68| for repulsion 
arising from going off-shell in the interaction is incorrect. Such repulsions do not exist. As 



demonstrated in Appendix F of [S4|, the physics does not depend on the way the kaon field 
interpolates. 

Schaffner et al j7(| point out the inconsistency of using the chiral vector and scalar 
interactions for the K~ mean field in the Walecka-type mean-field formalism for the nucleons 
(or nucleons and hyperons). The chiral Lagrangian gives V K + N = \V~nn rather than the 
\Vnn we discussed in eq. In Sakurai's work and in the schemes treating the vector 

mesons as gauge particles in an effective theory, which we discussed in Section 5, the uj- 
meson should couple to the baryon number lodged in nonstrange quarks. The p-meson 
couples to the isospin of the nucleon and this, plus the coupling to the pion, is responsible 
for the Weinberg- Tomozawa term in pion-nucleon scattering. The universal vector coupling 
can be obtained from the KSRF relation (23), and from it, the /oA^iV-coupling: 

^ = ^NN = . 70(5) (84) 

the (g/2) factor entering here because the isospin of the nucleon is 1/2. Via SU(3) symmetry, 
which is incorporated in the chiral Lagrangian, we obtain for the ujNN coupling 

= 9 ^nn ra 6 .34. (85) 

47T 47T 

This is substantially smaller than the loNN coupling constant employed in Walecka-type 
mean- field calculations: 

a 2 

S 10 ~ 12. (86) 

(The chiral uiNN coupling of (^) would give a vector mean field of only 178 MeV at nuclear 
matter density.) For the chiral calculation, we are justified in using the density dependence 
( |82D for /t,- derived from chiral considerations. Inserting ( |8"2| ) with the KSRF relation, we 
find 

9uNn(p = go) ~ g (87) 

4-7T 
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In other words, the connection between the chiral mean fields and those used in Walecka 
mean-field theory, which generally works at p ~ po, can be made if and only if BR scaling 
is taken into account. Neither theory employs form factors, so that comparison of coupling 
constants is straightforward. 

Calculations to date of kaon condensation in the Walecka mean field formalism have, 
indeed, been inconsistent in that they use empirical mean fields in baryonic interactions, 
but scalar and vector mean fields for the kaon from chiral Lagrangians without scaling /*. 

Our argument leading to ( p^ ) is difficult to employ for densities exceeding p . Even 
if the pion mass remains essentially unchanged, so that we can employ (^), the linear 
extrapolation (^) cannot be extended to higher densities. On the other hand the relativistic 
VUU calculations which employ Walecka mean fields work well up to several times nuclear 
matter density. 

We can turn matters around, and use the Walecka mean fields determined by the VUU 
calculations. But then the V^+n should be obtained as in (|79|), by the scaling indicated by 
quark counting, from Vnn- 

What to do with the Kaplan-Nelson scalar term is more intricate, since its orgin is in the 
chiral symmetry breaking part of the chiral Lagrangian. The factor F = (1 - 0.37uj 2 k /m 2 K ) 
which arises from virtual pair corrections at the same order in chiral counting as the Kaplan- 
Nelson term, must be included in an ad hoc fashion in the relativistic mean field formalism. 
Aside from this factor, scaling @ by the factor (0.6)- 1 of (||) gives S K n ~ -107 MeV 
and our determination of Skn from the lattice calculation of (iV|ss|iV"), as described in 
the Appendix, would give Skn ~ —120 MeV. Assuming quark scaling, this would give 
S n N ~ —360 MeV for the nucleon mean field, not far from that in Walecka theory. 

Although the scalar attraction on the kaon does not arise from a Walecka-type mean 
field mechanism, since the kaon is a Goldstone boson, use of the Walecka theory would give 
a rather similar effective scalar mean field. 

The factor F = (1 — 0.37uj k /m K ) is, however, important. For the if + -meson for which 
the ujk+ increases above mx for low densities at least, this factor cuts down the scalar 
attraction and makes the overall mean field repulsive, as is observed in i^ + -nucleus scat- 
tering. In the case of if - -mesons, this factor moves towards unity as ojk decreases with 
increasing density. Thus, it is not very important in determining the critical density for 
kaon condensation and the equation of state in the kaon-condensed phase. 

The importance of the above development is that it enables us to directly use information 
obtained from the Bevalac and SIS heavy-ion experiments, as described by tbe relativistic 
VUU transport in the calculation of kaon condensation. Thus we have empirical determi- 
nations up to p ~ 3po- As noted, from the subthreshold K + production experiments, we 
can already say that the vector and scalar mean fields, the latter corrected by the factor 
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(1 — oj 2 k+ / m? K ) , are roughly equal at p ~ (2i — 3) p$. 

Within the Walecka-type mean field description, we can now understand the large at- 
traction, —200 ± 20 MeV at p = 0.97/9o, for ET _ -mesons, found by Friedman, Gal and 
Batty J90|, This is just the sum of scalar and vector mean fields for the K~ . The number 
should not be surprising, since the same contribution of mean fields comes in the spin-orbit 
interaction for the nucleon in Walecka theory, and the above interaction for the K~ would 
imply that Snn + Vnn — 600 MeV which is big, though somewhat smaller than needed 
for agreement with experiments. (Of course we should not forget that our virtual pair 
correction cuts down the K~ interaction somewhat, ~ 10% at p ~ po-) 



6.3 "Cool" Kaons 



The 14.6 GeV 28 Si + Pb — > K + (K~~) +X preliminary data [79] show cool components, with 
effective temperature of 12 MeV for K + and 10 MeV for K~, which cannot be reproduced in 
the conventional scenarios employed in event generators. The latter give kaons of effective 
temperature ~ 150 MeV. It is clear that some cooling mechanism is necessary to produce 
the cool kaon component. It is also clear that the above vector repulsion must essentially 
be absent for the temperature relevant for these experiments, as we discuss. 

In an interesting article, V. Koch [8C] has shown that given the attraction of (^) - 



although his results are not sensitive to the precise amount of attraction as long as it is 
as large as ( jSP) ) -a cool kaon component can be reproduced. See Fig.[l4|. Aside from the 
attractive interaction, it is necessary that the fireball expand slowly. The slow expansion 
results because the pressure in the region for some distance above T x sr is very low Q, 
the energy in the system going into decondensing gluons rather than giving pressure. This 
results in an expansion velocity of v/c ~ 0.1. In the case of 15 GeV/N Si on Pb transitions, 
the fireball has been measured |^lj through Hanbury-Brown-Twiss correlations of the pions 
to increase from a transverse size of Rx(Si) = 2.5 fm to Rt = 6.7 fm, nearly a factor of 3, 
before pions freeze out. With an expansion velocity of v/c ~ 0.1, this means an expansion 
time of ~ 25 — 30 fm/c. (The full expansion time cannot be measured from the pions which 
occur as a short flash at the end.) Thus the fireball expansion lasts a very long time. 
As long as the expansion is slow, the adiabatic invariant 

I = J p • dr (88) 

remains a useful concept for the kaons. If the radius of transverse expansion r increases by 
a factor of 3, then the momentum will decrease by the same factor, and the energy by a 
factor of ~ 9. Thus, as long as the expansion is slow, the kaons will be greatly cooled, say, 
by the factor of ~ 9. 

Before proceeding further, we pause here to elaborate on the slow-expansion scenario 



for which the low pressure is essential. In Fig. 13 we give the pressure and energy densities, 
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Figure 13: Entropy, energy and pressure densities deduced by Koch and Brown j|] from the 
lattice gauge calculations of Kogut et al . 
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along with the entropy, deduced by Koch and Brown || from Kogut et al's lattice simula- 
tions § . Strictly speaking, the energy and pressure cannot be split up into a part from the 
quarks and a part from the gluons, for such a division is not gauge-invariant. So we should 
consider only the sum. 

What we see from Fig. pl| is that the pressure is only slightly above zero in the region 
where the entropy increases rapidly, i.e, in the region of the phase transition at T c ~ 
140 MeV. As noted in ||, this small pressure results because the energy chiefly goes into 
decondensing the quarks and gluons, thereby producing almost no pressure. This can be 
expressed in terms of an effective bag constant B for the transition. In order to see how 
this goes, consider the simplified case of a transition from pions to a quark/gluon plasma 



[38]. For simplicity, we neglect the pion pressure and energy, which can easily be corrected 



for at the end. The energy-density of a quark/gluon plasma is 

e QG = ^ 7f2T4 + jBe //' 

Pqg = ^ 7f2T4 - B eff- 

Now the phase transition occurs as soon as the Pqg can be brought positive. From lattice 
calculations, we know that T c 140 MeV. From Pqg = at T c we find 

B X J f ) = 199MeV, (90) 

midway between B x $r of eq.(||) and B g i ue of eq.(|5|). Koch and Brown show that this 
"half-way house" results because only ~ | of the gluon condensate is "melted" across the 
transition region. Note that if we were to calculate T c from first principles, we would need 
to know this latter fact. In other words, the determination of -B e //> that fraction of B which 
is melted, must first be found. By using the lattice results, we have turned the problem 



around, obtaining B e ff. Equation (89) makes it clear why the pressure is so small in the 
region of T c . 

The smooth transition that we are finding implies that instead of having the "melting" 
occurring at a fixed T = T c , it occurs over a region, determined as AT ~ 20 MeV by Koch 
and Brown |J, although this number may change as nonperturbative effects are included 
in the f3 function. 

Estimates of the energy densities reached in the Brookhaven AGS collisions are quite 
model-dependent, but generally indicate that the highest temperatures are not more than 
~ 170 MeV, and probably are less than this. Consequently, the systems do not reach 
temperatures much above T c . 

We should remark that while the relationship between hadron masses and the quark 
condensate, investigated in || does not involve knowing the temperature, the width of the 
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phase transition (~ 20 MeV) does depend upon the asymptotic scaling relation being effec- 
tively valid. One can approximately check the validity of this relation as follows: Asymptotic 
scaling can be expressed as 

8tt^ 2 k 

H_ jV F in- 91 

T 3 Aqod 

giving the relation between the color gauge coupling g and the momentum k. For a lattice 
of 6 x 12 3 , Kogut et al |j| found the transition at [3 C = 6/ ' g 2 c 5.34, whereas for 8 x 16 3 , 
Gottlieb et al p^] find j3 c = 5.54. Assuming the scale to depend chiefly on the number of 
time slices, one can get rid of Aqcd using the two simulations, in order to express 

^ = ^(ll-|nj,)hi(|)«.21. (92) 

This "theoretical" 50 c , which assumes asymptotic scaling, is to be compared with the 8j3 c ~ 
0.2 found in the lattice simulations. The good agreement indicates that the asymptotic 
scaling relation should be adequate for determining temperature differences. 

Now in expansion of the system following maximum temperature, the expansion will 
"stick" for a long time (estimated to be 25 - 30 fm/c) just in the region of T c , where 
the pressure is nearly zero. Entropy density is decreased not by expansion, but by the 
hadrons, which are formed in this region, going back on-shell and the number of heavy 
hadrons then decreasing because of small Boltzmann factor e~ mH / T . This picture seems 
to be supported by the work of Shuryak and Xiong |84j who find evidence in the excess 
photons and dileptons in the SPS collisions (200 GeV/nucleon collision at CERN) for a long 
"mixed" phase of r ~ 30 - 40 fm/c. 

Now returning to the problem of the cool kaons, the important point to note is that the 
kaons move through chirally restored matter, in which the hadronic (as opposed to gluonic) 
vector interactions are small or zero during the 25 — 30 fm/c period. Kaons which do collide 
strongly inelastically will be removed from the cool kaon component and join the thermal 
kaons. 

Given the cool kaon results, the collaboration E814 - now E877 - carried out a study 
of this component in Au+Au collisions at 10.8 A GeV/c at the AGS|53|]. Both spectra 
exhibited a cool, nonthermal component but with apparent temperatures in the ranges of 
~ 50 - 100 MeV (depending on the rapidity), substantially higher than ~ 10 - 12 MeV 
found for Si+Au. This higher effective temperature was predicted by Koch|8(J. The Au 
system is larger, and the kaons have more chance to hit something and be scattered out of 
the cool component. What was unexpected was that, with decreasing Pt, the spectra began 
to drop in the case of the K + spectrum but not for the K~ one. (Following the presentation 
of the work[^] at the Quark Matter '95 meeting, C.Y. Wong suggested that the drop in 
K + spectrum with P t decreasing below 140 MeV was due to the Coulomb interaction which 
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Figure 14: Kaon spectrum from a full transport calculation for 14. 6A GeV/c Si + Pb collisions [j79[. 
The full line is the result including mean fields for baryons and kaons, the dotted line the result 
for kaon mean field and no mean field for the baryons and the dashed line the pure cascade result. 
The calculated results are in arbitrary units. The nucleoli vector and scalar couplings of the mean 
field were taken to be gy — 5.5 and gs = 9.27, the vector coupling being about half of the Walecka 
mean- field model, so that the kaons - which felt 1/3 of the nucleon mean fields - experience an 
attractive potential of Uq ~ 50 MeV. The field energy |ms</>|, where <ps is the scalar field, plays the 
role of the bag constant (^). Its value, taken to be about half of the B of (||), was adjusted so that 
the pressure was low in the region of the phase transition, in accordance with the lattice results of 
Fig|| 
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would be substantial in the large Au system. Work by Koch, in a schematic model, showed 
this to be a plausible explanation. Further work with full transport by Koch is in progress.) 

It can be seen that the discovery of the cool kaon component tells us that the main part 
of the system is chirally restored (in the form of the Georgi vector limit), for a long time of 
~ 25 — 30 fm/c, until the kaons freeze out. They freeze out well ahead of the pions, which 
equilibrate with the nucleons down to a density, of pf ~ 0.38/JO) because of their weak 
interaction. 

Since masses of particles go back on-shell in free space before the particles reach the 
detector, the bulk of experiments is insensitive to what the masses were inside the hot and 
dense medium. Only carefully designed experiments will tell us that, and the cool kaon 
component is the most important of these to date. 

7 Discussion of Chiral Restoration with Temperature 

We have seen in the last section that the behavior of the quark number susceptibility 
with temperature, as calculated in lattice gauge calculations, shows that the hadronic vector 
coupling disappears as T moves upwards through T x $r, and that the perturbative color 
gluon exchange describes the susceptibility well above T x sr, as argued by Prakash and 
Zahed |52| . Somewhat surprising is the fact that the perturbative description, which gives 
a 1/T 2 behavior in the difference x(T) — x(0) between the susceptibility and that for free 
quarks, sets in just above T x sr; i.e, for this purpose, asymptotia is T j> T x sr. Note that 
the screening mass of the /9-meson goes to 2ttT, its asymptotic value, as soon as T reaches 
T x sr- 

As noted preceding eq. (pfl) , the hidden local gauge coupling g, at one loop order, scales 
to zero as g(/i) ~ (In /x) _1 , i.e, the Georgi vector limit. As shown by Harada and Yamawaki 
pi] , there is an ultraviolet fixed point tz = 0. Of course, the one-loop calculations are most 
certainly unreliable for learning what happens near the phase transition, where the higher 
loop effects and nonperturbative effects presumably pile up. Even so, the qualitative picture 
suggested by this one-loop calculation appears to be correct. 

It may be that the increase in /* = /* as T goes through T x sr found in the last section 
is an artifact of the spacelike propagated wave function. With timelike propagation, which 
has not been up to now carried out on the lattice, /* and /* may go to zero as T — > T x sr, 
as suggested by the BR scaling Q. What one does learn from the spacelike propagated 
wave function is that the quarks in the p-meson behave, to lowest order, as free quarks, 
each with thermal energy of 7rT, with corrections of ~ a s which give the Bethe-Salpeter 
amplitude for the p. It is, in fact, this basically free quark kinematics which leads to 
the equality /* = /* for the spacelike wave functions. Wave functions propagated in the 
time direction will, presumably, have this same basically free quark kinematics. Indeed, 
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in four-flavor lattice calculations of the p-meson [41], the temporal correlators above T x $r 
were found to be consistent with a pair of quarks that are more or less free. Modulo the 
caveat with four-flavor lattice calculations mentioned before, given free quark kinematics, 
one would expect /* = /*, even if /* — > with chiral restoration. In the possible case of 
/* — > as T — > T x sr, the Georgi limit could be realized either just before or simultaneously 
with the /* going to zero. Just what happens in the phase transition may be complicated 
to unravel. Nonetheless, the lattice calculations of the susceptibility do show that, to the 
accuracy we can analyze them, g — ► as T — > TySR, indicative of the Georgi vector limit. 

While other viable explanations may be found in the future, as far we know, the soft 
kaons found in the E-814 experiment can be explained only if the vector-meson coupling is 
essentially absent at the relevant temperature, and this is a strong support for the Georgi 
picture in which the hidden gauge coupling constant "melts" at the temperature T x $r. 

Since Sakurai [82] we have known that vector dominance and universal vector coupling 
worked well in describing vector interactions. It is in endowing the vector mesons an induced 
gauge structure [jn| that a contact (albeit an indirect one) with QCD is made. As pointed 
out by Georgi [17], the hidden local symmetry is useful because it keeps track of powers 
of m p /A x , where A x is the scale of the effective chiral theory, in the situation where the 
vector-meson mass m p is in some sense small. Indeed the successful KSRF relation follows 



simply [17] if the vector mesons as well as the pseudoscalars are both light, so that we can 
stop with a Lagrangian with the lowest (that is, two) power of the derivative dictated by 
chiral invariance. This suggests rather strongly that the notion that the vector-meson mass 
becomes small as T — > T x sr is natural within the framework of effective chiral field theory. 
Remarkably, this seems to be supported by lattice gauge calculations at T ~ T x sr. Now the 
gauge symmetry in the hidden symmetry scheme is an "induced" gauge symmetry lodged 
in hadronic variables. In this sector, the fundamental color gauge symmetry is not visible. 
It is the induced flavor one that is seen. What we observe is then that as T goes towards 
T x sr, the induced gauge symmetry gives way to the fundamental gauge symmetry. What is 
surprising is that this changeover seems to take place suddenly, with increasing temperature, 
the effective hadronic gauge symmetry applying for T < T x $r, and the fundamental color 
gauge symmetry being realized perturbatively for T > T x $r. 

We have suggested that the Georgi vector symmetry with k = and g ^ is also 
relevant for nuclear physics. It is the existence of such a symmetry that allows one to 
linearize the non-linear chiral Lagrangian of current algebra to the linear sigma model 
which justifies in some sense the separation of the scalar x held, the low-frequency part of 
which giving rise to the previously obtained medium-scaling of effective chiral Lagrangians 
of ref.pj. Going beyond the tree approximation with this scaled chiral Lagrangian in hot 
and/or dense matter is an open problem which has attracted little attention up to date. 
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While consistent generally with both observations and theoretical prejudices, our discussion 
cannot be considered solid until higher-order calculations can be systematically carried out 
and compared with experiments. A small progress made in this direction is discussed in 
ref.|§. 

Finally we should mention that an investigation of effective hadronic interactions medi- 



ated by "instanton molecules" carried out by Schafer et al [33| provides support to our thesis 
that chiral phase transition involves Georgi's vector limit. In this work, it is found that 
all coupling constants in a NJL-type effective Lagrangian can be specified for T J> T x sr in 
terms of a single parameter G, signaling a swelling symmetry. At temperature T > T x $r, 
the instantons can be taken to be completely polarized. In this case, the interaction in the 
longitudinal vector-meson channel becomes equally strong as the attraction in the scalar- 
pseudoscalar channel. Transversely polarized vector mesons are found to have no interac- 
tion. In detail how this comes about is as follows: The instanton molecules survive the chiral 
symmetry restoring transition; they leave chiral symmetry unbroken. According to Koch 
and Brown M, about half of the gluon condensate remains for temperature T > T x sr. This 
is assumed to reside in the molecules p3| . Quarks coupling through the instanton molecules 



experience, for randomly oriented molecules, an interaction 

2C 1 

6H = J^{(^^(x)) 2 + -(^r^(x)) 2 + • • •} (93) 

where we show only the terms that act in the pionic and p-meson sector. The \ in front 
of the vector interaction comes from averaging the molecules over random directions. In 
this case, it is seen that the interaction in the p-meson channel is only | as large as pionic 
channel. If the instanton is completely polarized in the time direction - and the authors of 
|83| give reasons why this lowers the energy of the system - the interaction in the p-meson 
channel changes to 

-^^{x)) 2 -» {i> l0 T^{x)f (94) 

since directions are no longer averaged over. The \ factor no longer appears, because there 
is no averaging over directions. This sets interaction in the pionic channel and that for the 
time-like p meson equal, decoupling the transverse p mesons. The time-like vector mesons 
are longitudinal in the sense that 

<f Pli = (95) 

where is the four-momentum of the p meson. We thus have a "realistic" dynamical 
model, giving the effective interaction in hadronic variables, which reproduces the Georgi 
vector limit under the assumption of a dilute gas of fully polarized instanton molecules and, 
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of course, the p meson going massless. p^l (In this model, the u and the p can be put into 
flavor U(2) as well.) Note that the above considerations apply to the region of T just above 
T c , where the colored gluon exchange is present. 



8 Chiral Restoration in Dense Matter in Stellar Collapse: 
Transition to Kaon Condensation and then to Quark Mat- 
ter 

In Section 6 we discussed effects of the attractive kaon-nucleon interaction, namely, the 
Kaplan-Nelson term, on subthreshold production of mesons in heavy-ion collisions of 
Au on Au, with 1 GeV/nucleon energy. These involve chiefly finite density effects, because 
the temperature is low, ~ 75 MeV, so that T is well below T c . We do not expect much 
effect of dropping masses, etc., because of the temperature. 

Our consideration of chiral restoration with high temperature could be checked by lattice 
gauge simulations; more correctly, many of our ideas about the restoration stemmed from 
lattice gauge results. Unfortunately, lattice calculations cannot be carried out for finite 
density - at least not to date - because in changing to imaginary time (Euclidean space- 
time) the baryon chemical potential becomes imaginary and the probability is no longer 
positive, so that Monte Carlo methods no longer work, at least not effectively in their 
present formulation. 

Effects from finite density in strong interactions have been reviewed extensively by 
Adami and Brown || and so we will not repeat this review here. A very interesting situation 
is however provided by the collapse of stars, in which there is sufficient time for strangeness 
conservation to be violated. We shall review recent developments on this here. 

For concreteness, let us sketch the scenario for Supernova 1987A. According to the 
scenario we shall pursue here, the core of the 18 M & progenitor collapsed, initially forming 
a neutron star, and blew off the outer material, in a time of ~ 4 seconds following collapse 
fl8ll . During this time the neutron star was stabilized by the trapped neutrinos. We know 
from the Kamiokande neutrino detector that the neutrinos came off for ~ 12 seconds; by this 
time they had dropped in energy below the ~ 8 MeV threshold for detectability. According 



to Burrows and Lattimer [87] only about half of the thermal energy is carried away by the 
neutrinos during this ~ 12 seconds. About half of the thermal energy remains and creates 
substantial pressure, which helps to stabilize the neutron star. The central density of the 
neutron star grows to several times nuclear matter density. 



* 22 Schafer et al |g3[ find that the molecules are not completely polarized (the polarization being ~ 70%), 
so this model may be taken only as a caricature theory which, in a limit not far away from the physical 
world, reproduces the results of the Georgi vector limit. 
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Once the electron neutrinos have left, following the explosion, the electrons in the core 
turn into if "-mesons and neutrinos, as we now explain. The neutrinos leave the core, again 
taking several seconds to do so. This cannot happen while the first set of neutrinos are 
trapped, because this would require 

W/f = HK = He~ Hu, (96) 

where the /x's are the chemical potentials, and as we shall see, it would be difficult to bring 
the kaon energy lok down this far. But when the \i v — > 0, it then becomes possible to bring 
•jJk down to the electron chemical potential \x e . 

Prakash, Ainsworth and Lattimer ]8S| give simple parametrizations of the neutron rich 
equation of state for the situation after the neutrinos have left. The electron chemical 
potential is then established by the thermal equilibrium condition 

e~ + p <-> v + n, (97) 

with the neutrinos more or less freely leaving the star, so 

(l e = H n - Hp- (98) 

The chemical potentials of the neutron and proton must be calculated, and \i e determined 
from these and the condition of charge neutrality. Although He depends somewhat on the 
compression modulus, etc., a good standard value for p e is 

He = H7MeF for p = p . (99) 

The He increases slightly faster than as p 1//3 with density, reaching > 200 MeV at about 3po 
for equations of state with a nuclear matter compression modulus of Kq = 200 MeV. These 
numbers schematize the results of Thorsson, Prakash and Lattimer |39| and of p3j. 

Now the energy of a If "-meson has been found to drop in nuclear matter. According 



to Friedman et al |K|, the real part of the -fT - -nucleus potential is —200 ± 20 MeV in the 
center of the 56 Ni nucleus at p = .97 pq. The K~ feels double the attraction from a proton 
as it does from a neutron through the vector interaction (|79|) and the same attraction for 



neutron and proton through the scalar interaction (8C). Thus, for neutron-rich matter, one 
can say that the K~ energy is decreased ~ 150 MeV at nuclear matter density po- What 
can happen is schematically described in Fig. 15. 

At the point where the lok line crosses the He h ne > it becomes energetically possible for 
the "electrons" to turn into K~'s accompanied by neutrinos 

e~ -► K~ + v e (100) 

with the neutrinos leaving the star. There will therefore be a second burst of electron 
neutrinos which should be detectable. Details (energy, time of emission, etc.) are being 
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Figure 15: Behavior of the K energy and of the electron chemical potential fi e as function 
of energy. 

worked out by A. Burrows, J.M. Lattimer, M. Prakash and collaborators in Stony Brook 
at the present time. 

We note that the large binding energies of Friedman, Gal and Batty fit in very well 
with the scenario of Brown-Rho scaling, with / — > /* in ( fTS|) and (0C| ) and as shown in 



Figs. 11-13 of the detailed calculations in [64]. One seems to need as much attraction as 
the theory can provide. 

At densities above the kaon condensation threshold, the K~ mesons of zero momentum 
form a Bose condensate, as discussed in [|89], H]. Sufficient binding energy is gained so 
that the equation of state is substantially softened, and the maximum neutron star mass is 



~ 1.5Mq as discussed by Brown and Bethe [91]. We will return to this matter at the end 
of this section. 

It has been recently shown by Brown and Weingartner [p!| that in the case of Supernova 
87A if a neutron star were present, then we would see it with a luminosity L of about 10 4 
times Lq, whereas the bolometric luminosity actually is L ~ 1O 2 L and can be explained 
as arising from radioactivity in the expanding shells. Thus their conclusion is that as in 
PH the core has gone into a black hole. Furthermore Bethe and Brown [93] show that from 
the fact that 0.05 M of Fe was produced, together with the pre-supernova evolution, it is 
possible to determine the mass of the compact object to be in the range (1.44 - 1.56) Mq, 
confirming the suggestion by Brown and Bethe f9l|| . 

Recently Keil and Janka | ]94| | have shown, based on the Glendenning mean- field equations 
p5|, that hyperons enter into the equation of state at densities p ~ 2pQ. They find a delayed 



explosion scenario similar to that of Brown and Bethe [91] for cores in the range 



1.58M < Mg rav < 1.72M©. (101) 
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Here we have used the binding energy of 

E = O.O84M (M/M Q ) 2 (102) 



of Lattimer and Yahil [96] in order to connect baryon number masses given by Keil and 
Janka with gravitational masses. Undoubtedly parameters could be adjusted to bring the 
lower limit down to the 1.5 M Q of Bethe and Brown ]^3|, which is required for observation. 
More extensive calculations with inclusion of exchange (Fock) terms p7| bring the maxi- 
mum neutron star mass above 2 Mq, the Fock terms significantly stiffening the EOS once 
hypersons are present. 

Introduction of hyperons into the EOS in Hartree approximation causes the electron 
chemical potential to saturate |M| at p e ~ 200 MeV at p ~ 2pQ. If hyperons enter at 
densities before kaon condensation, they may hinder the condensation by bringing p, e down. 
On the other hand, they sufficiently soften the EOS that the system moves rapidly to 
a higher density, which brings the K~ energy ujk down, and these effects are likely to 
compensate. As noted, linear extrapolation of the binding energy of a K~ in nuclei to higher 



densities, or inclusion of Brown- Rho scaling 64], gives p c for kaon condensation only slightly 



above 2p$. Once finite temperature is included, a substantial number of thermal K -mesons 
are present, so the distinction between introduction of hyperons and kaon condensation is 



likely to be smeared out. Recently Ellis, Knorren and Prakash[98] have introduced hyperons 
into the dense matter, largely confirming the calculations of Glendenning[95] and of Keil and 
Janka Ellis et al find that the introduction of E~'s brings down the electron chemical 
potential p e and delays the kaon condensation until higher densities. We believe this to be 
incorrect for the following reasons: (1) E~'s do not experience the same attraction in nuclei 
as the A's do|}9), although this is assumed in the quoted calculation; a possible explanation 



is that QCD sum-rule calculations [10C] find substantially more repulsive interaction on the 
£~ than on the A, in contrast to quark scaling which predicts these to be the same; (2) 
the Ellis et al calculations use Walecka mean fields for the hyperons, but interactions from 
chiral Lagrangians for the K~ mean fields. As noted at the end of Subsection 6.2, this is 
inconsistent. Given the Walecka mean fields for hyperons, the K~ mean fields should be 
a factor ~ 1.6 larger than employed. This not only makes kaon condensation competitive 
with the introduction of hyperons, but also means that it will probably precede any such 
introduction. To put it more directly, we believe that kaon condensation will win out over 
introduction of hyperons. Even if the latter did indeed enter first at lower densities, it would 
soften the EOS sufficiently in order to lead quickly to kaon condensation. 

Aside from substantially softening the EOS, so that a compact object of mass 1.5Mq can 
go into a black hole, the strangeness condensed EOS has an important consequence, shown 
in Fig.[ji] [101]. While the neutrinos are trapped, the leptons (neutrinos and electrons) 



produce substantial pressure, so that the maximum mass given by the solid line in FigJT^ is 
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Figure 16: The solid line shows the short-time maximum neutron star masses as function 
of central density, with an assumed lepton fraction of Y\ = 0.4, consisting of electrons and 
trapped neutrinos. The dashed line shows the long-time maximum cold mass, after the 
neutrinos have left. Note that the highest n cen t is ~ 12po- 



stabilized. Once the neutrinos leave, the electron fraction Y\ drops to a low value, and since 
most of the lepton pressure has been carried away by neutrinos, the maximum mass that 



can be stabilized is substantially less, 1.5M in the case of Fig. 16. Consequently, compact 
objects with masses in the range 

1.5M < M < 1.7M (103) 

will be stable during the time of neutrino emission; i.e, long enough to explode, and will 
then go into a black hole. Brown and Bethe [91] estimate that this delayed drop, after 



explosion and return of matter to the galaxy, into a black hole will take place for most stars 
with main sequence masses 

18M < M < 3OM . (104) 

Supernova 87A lies near the lower end; in fact, it essentially sets the lower end. In 
PSR 1913 + 16 the pulsar, the larger of the neutron stars in the binary, has mass 1.44M . 



Evolutionary calculations [102, 102] find that a 1.44M neutron star results from a 5 - 6 
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Figure 17: A very soft EOS with Kq = 130 MeV so that the maximum neutron star mass 
for cold matter (dashed line) is 1.5M@ as in Fig.|l~6|. The short time maximum mass, given 
by the solid line, is now lower, even though the neutrinos are trapped with Y\ = 0.4 
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Mq helium star, corresponding to a main sequence mass of 16 - 18 M & for the progenitor. 
Thus, the pulsar in 1913 + 16 is very close to the maximum mass for a neutron star. 

Not only do we believe that the Supernova 87A ended up as a black hole, but this was 
also likely for CAS A, a supernova explosion that took place in the 17th century. From the 
abundances of O, Mg, Ne in knots in the supernova remnant, the masses of the progenitor 
can be deduced to be M ~ 20M & . No compact object is found in the center of the remnant, 
although searches, with great sensitivity, have been carried out. 

We shall now explain why the delayed explosion does not result from the standard neu- 



tron star scenario. This situation is shown in Fig. 17. In this case the short-time maximum 
mass, with neutrinos trapped, and with Y\ = 0.4, is lower than the cold mass, given by the 
dashed line, after the neutrinos leave. The point is that, in the standard scenario, as the 
neutrinos leave, the original nuclear matter is converted to neutron matter, and neutron 
matter is much "stiffer" (that is, the pressure is higher) than nuclear matter. Thus, even 
though the pressure is lowered by the neutrinos departing, it is increased somewhat more, 
by the protons changing to neutrons. Thus, if the compact object is initially stable, it will 
remain stable, and will not go into a black hole. 

Previous to the scenario involving the kaon condensed EOS, it was suggested that stars 
heavier than ~ 25M may leave black holes [104], and also that such stars might first 



explode, exhibiting light curves of Type II supernovae, and then collapse into black holes 
|105| ], The compact core was, for certain range of masses, to be stabilized by the thermal 
pressure during the period of Klein-Helmholtz contraction, long enough to carry out nu- 
cleosynthesis, then going into a black hole after cooling and deleptonization. In terms of 



Fig. 17, this scenario meant that the thermal pressure had to be sufficient, so that when 
thermal effects were added to the (short-time) solid line, it came close to the dashed line. 
Then, the compact object could have a mass lying above the dashed line (and below the 
solid line) and be stable for some time, before it went below the dashed line as it cooled. 

Table 1 

Energy gain AE in MeV, chemical potential px hi MeV, proton fraction x and 
electron fraction x e as function of the density u = p/po for a kaon condensed EOS 
with u r = 4.2. 



u 


AE(MeV) 


PK(MeV) 


X 


x e 


4.2 





256 


0.20 


0.11 


5.2 


-10 


199 


0.34 


0.03 


6.2 


-35 


142 


0.43 


0.01 


7.2 


-71 


94 


0.48 





8.2 


-112 


55 


0.50 
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The trouble with this scenario, of thermal pressure stabilizing the compact object, is 
that the thermal pressure can stabilize only a small additional mass. Detailed calculations 
by Bombaci et al [106] give only an additional ~ 2-3% increase in stabilized mass, less than 
the distance between the solid and dashed lines in Fig. 17. Similar results are obtained by 



Keil and Janka [94] with their delayed explosion scenario in the EOS including hyperons. 
Thus, it is clear that in most cases, thermal pressure will not stabilize the compact object 
for some time, with the object later going into a black hole; rather, it will remain stable if 
it is so initially. 

We see that a lot of observational evidence is explained with our kaon condensation EOS. 
Note that in this scenario, stars go "strange"; i.e, acquire a lot of strangeness already in the 
hadron sector. They do not, at the densities of 2 - 4 po, go to strange quark matter. The 
transition from kaon condensation to strange quark matter can, however, be constructed 
[107]. This was made under quite conservative assumptions of a small Xktv = l-3m T and no 
Brown-Rho scaling. We give in Table 1 the results of the calculation by Vesteinn Thorsson 
in his 1992 Stony Brook thesis. With small differences, these are the same numbers as in 
Table 3 of ref.§l|. 



Of chief concern to us is the softening of the EOS by kaon condensation. In Fig. 18 



we plot the baryon number chemical potential vs. pressure. This is a convenient plot 
for investigating the transition to quark matter. In the following we use the procedure of 
Bethe et al gof]. The PAL curve is taken from Prakash et al |88|] . It has a compression 
modulus Kq = 180 MeV, and the potential energy part of the symmetry energy rises 
linearly with u = p/po- The PAL21 curve from 88] has a maximum neutron star mass of 
M max = 1.72Mq. The curve including the kaon condensation with characteristics shown 
in Table 1 is plotted in the lower dashed line in Fig. 18. In addition, the results for quark 



matter [ 108 ] are given, but with strange quark mass m s = 200 MeV included. (In flOSj l, m s 
was set equal to zero.) There a s was scaled as 

a s (k F ) = 2.2kP/k F (105) 

where kp is the Fermi momentum at nuclear matter density, kp the Fermi momentum at 
the density considered. This gave a rapid decrease in coupling constant, which Bethe et al 



[108] considered appropriate for the nonperturbative sector. 

It should be noted that the M.I.T. a s of 2.2 at p = po, even decreasing as rapidly as ( 105 ), 
gave an EOS which lies well above the PAL21 curve, so there is no hope of joining PAL21 to 
it, especially if PAL21 is decreased even further by kaon condensation. (A "conventional" 
compression modulus of Kq = 210 ± 30 MeV [74] would give an only slightly higher curve 
than PAL21.) 

Quark EOS's for a s = 1.1 and 0.55 are also shown. The Pq + B is plotted for the 
quark/gluon phase, so that introduction of a bag constant (which does not affect the chem- 
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Figure 18: Curves of chemical potential vs. pressure for hadronic matter (dashed curves). The 



upper dashed curve representing PAL21|88 has Kq = 180 MeV, a symmetry energy for which the 
potential energy rises linearly with density. The lower dashed curve results when kaon condensation 
is included. Solid lines for quark matter are plotted for pq + B, instead of the pressure, but B will 
be taken to be zero. The upper solid line is for a s — 1.1; the lower one, for a s — 0.55. On each line, 
the lowest black circle or square marks the density 2po, and each successive dot or square indicates 
a density higher by po. 
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ical potential p) can be made by shifting the Pq + B curve to the left by the amount B, 
so as to obtain Pq. As outlined in [||, we believe the bag constant, the B x sb of eq.@, 
to go to zero at chiral restoration. Note that the quark/gluon curve for a s = 1.1 is easily 
made equal to PAL21 over a wide range of pressure about Pjv ~ AOOMeV/ fm 3 in this way. 
Similarly, the quark/gluon curve with a s = 0.55 can be made equal to the PAL21 with kaon 
condensate curve, in the region of pressures which correspond to densities of ~ 8-10po for 
hadronic matter. (The density is not continuous at the transition.) We see immediately 
that with the M.I.T. value for a s used by Bethe et al [105], there is no hope of making 
a transition to quark matter; the quark matter EOS has a much higher energy than the 
hadronic one, especially when kaon condensation is introduced into the latter. 



Fahri and Jaffe [10£] have taken the position that the a s for the dense matter properties 
may be substantially smaller than 2.2, and have investigated the range of a s used in the 
two curves of our figure. Given these small a s , it is possible to discuss the transition to 
quark matter. 

With our conservative choice of the kaon condensed EOS, without Brown-Rho scaling, 

fc (o) 

and the a s of 0.55-^, a smooth join to quark matter can be made at a density of ~ 8- 
10po- Whereas a s may not drop as rapidly as kp /}cf, an a s of ~ 0.25 at p ~ 10po is not 
unreasonable, although we have no way at present of determining a quantitative value at 
such high densities. We can thus say that, with a rather conservative kaon condensed EOS 
(with p c ~ 4po)j there is the possibility of a smooth cross-over transition to quark matter 
at the upper end of the densities obtained in compact objects, say, p ~ 10po- With a softer 
EOS, obtained with Brown-Rho scaling, with kaon condensation taking place at p <i 3po, 
the transition density to quark matter would be substantially higher. 

It is daring - and perhaps foolhardy - to extrapolate to densities ~ 10po, but we believe 
the replacement of electrons by X~-mesons at high densities to be a new idea with a solid 
foundation, which qualitatively changes the conceptual situation in compact star matter. 
We may expect the transition to quark matter at high densities to take place similarly to 
the transition with temperature to quark/gluon plasma, as discussed in Section 4; namely, 
the transition will be smooth and gradual. At the lower densities it will be more convenient 
- and cleverer - to use hadronic variables, but at the higher densities the quark language 
will become more efficient. As argued, we may expect the light-quark vector-meson masses 
to go to zero, with increasing density, and chiral symmetry to be realized in the Georgi 
vector limit preceding chiral restoration. 

Given the kaon condensation phase transition, our EOS is so soft that it is difficult to 
stabilize stars of known masses. Indeed with the PAL21 EOS including kaon condensation, 
M max = 1A2Mq, although raising Kq to the more conventional value of 200 MeV allows us 
to obtain M max = 1.5-M©. A slightly higher Kq will be needed if the kaon condensed EOS 
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with Brown-Rho scaling is used, since this EOS is substantially softened by this. 

Observed neutron star masses are shown in FigJlfl. Precisely the lower limit of the 



measured mass of Vela X-l lay below 1.5-Mq. New observation by Van Kerkwijk et al [ 111 ] 
found that the observed velocities in Vela X-l deviate substantially from the smooth radius- 
velocity curve expected from pure Keplerian motion. The deviations seem to be correlated 
with each other within one night, but not from one night to another. The excursions suggest 
something like pulsational coupling to the radial motion, and make it difficult to obtain an 
accurate mass measurement. The lower limit for the mass of the compact object in Vela X-l 
is now found to be 1.43M at 95% confidence level, or 1.37Mq at this confidence interval 
around the most probable value. Consequently, Vela X-l is no longer a big problem for our 
M max = 1.5M . 

It is striking that well measured neutron star masses lie below 1.5Mq. (See Fig.|l9|,) 
However, the central value of the compact object in 4U 1700-37 lies at 1.8M@, although the 



error bars encompass 1.5Mq. Brown et al [103] give arguments that this compact object 



could be a low-mass black hole. (The principal aim of this paper is to show that stars in 
binaries can have much larger masses than the main sequence mass of 18M Q of 1987A, 
and still end up as neutron stars. In other words, in binaries, because of the specifics 
of mass transfer, stars can evolve in quite a different way than single stars evolve. Since 
determination of masses has generally been carried out in binaries, this explains why it 
was not earlier recognized that a star with main sequence mass as light as ~ 18M (the 
progenitor of 1967A) could go into a black hole, and it was - and still is - considerable 
surprise that it probably did.) 

In the past it has generally been thought that the reason accurately measured neutron 
star masses lie at 1.44M and below is evolutionary in nature. Large stars collapse when the 
iron core exceeds the Chandrasekhar limit, ~ 1.25— 1.5Mq, depending on the main sequence 
mass of the star. In the past literature, accretion has been assumed to proceed only up to 
the Eddington limit 

M Edd = 1.5 x 10- s M Q /yr (106) 

and there are relatively few situations where the compact core would be expected to accrete 
more than ~ 0.1M@ at this rate. (Even in the very old millisecond pulsars, only ~ 0.1- 
0.2Mq is estimated to have been accreted.) However, Chevalier ]112| ] and Brown [ |113| ] have 
shown that in the common envelope phase of binary pulsar evolution, accretion can proceed 
at hypercritical rates 

M > lQ A M Edd . (107) 

Thus, if evolutionary history determines their mass, neutron stars of higher mass than 
1.5M Q should exist, but none have been so far observed. Therefore, there should be an 
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Figure 19: Measured masses of 17 neutron stars from Arzoumanian et al [11C], with the 
lower limit on the mass of Vela X-l from Van Kerkwijk et al 1 1 1 1 ] . Objects in high mass 
X-ray binaries are at the top, radio pulsars and their companions at the bottom. 
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intrinsic limit on the mass of a neutron star, such as the one we find. It is difficult to 
explain the existence of such an intrinsic limit without softening of the EOS through a 
phase transition and we are proposing that kaon condensation is the key mechanism for it. 
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8.1 Appendix: The Kaplan-Nelson Attraction 



The attractive scalar mean field potential (p0[), used first by Kaplan and Nelson ]60[] for 
kaon condensation, depends linearly on the KN sigma term Y^kn, which in turn depends 
on the strangeness content of the nucleon (iV|ss|iV). This is usually parametrized by 

= 2(N\ss\N) 
V (N\{uu + dd)\N)' { ' 



In the most extensive lattice gauge calculation to date, Liu finds [114] 

y = 0.33 ±0.09. (109) 

Given this, we can estimate S^at using the information on the ttN sigma term E^jy- Since 

Zrn _ (m s + m u )(N\(ss + uu)\N)) 

X nN ~ (m u + m d )((N\{uu + dd)\N)) , [ ' 

taking the value for m s of 

2m s /{m u + m d ) « 29 (111) 

from Bijnens et al [115], y = 0.33 and T, n ^ = 45 MeV, we find 

Y, KN ^ 450 ±30 MeV (112) 

as the best current estimate of T>kn- This is slightly larger than the 2.83 m n found by Lee 
et al [64] from fitting the KN scattering amplitudes. We believe that the wild fluctuations 
in Tjkn used in the literature have now settled down. 
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Using Liu's value |TT^| of 

(N\uu + dd\N) = 8.22 ±1.1, (113) 

we find as central value 

(N\ss\N) = 1.36. (114) 

With m s = 174 MeV, this would mean that m s {N\ss\N), the contribution to the nucleon 
mass from the explicit chiral symmetry breaking in the strange sector, is 237 MeV which is 
sizable. 

In chiral Lagrangians, the explicit chiral symmetry breaking in the strange sector is 



parametrized by a coefficient denoted a 3 in [ 116 1 . Politzer and Wise employed rather dif- 



ferent values, a^m s = 310 MeV corresponding to a large (iV|ss|iV) and a^m s = 140 MeV 
corresponding to a small (N\ss\N)^^_. Using our above central values for m s and (N\ss\N), 



we find, using the relations given by Politzer and Wise, 

a 3 m s ps 245 MeV (115) 
with an estimated uncertainty of ~ 10%. 



Note that the T,kn = 450 MeV of eq.( 112 ) would give a coefficient of 72 MeV, rather 



than 64 MeV, in eq.(pO|). This difference is not significant, given the listed uncertainties. 
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